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We derive underdamped evolution equations for the order-parameter (OP) strains of a ferroelastic material 
undergoing a structural transition, using Lagrangian variations with Rayleigh dissipation, and a free energy as a 
polynomial expansion in the N = n + N op symmetry-adapted strains. The N op strain equations are structurally 
similar in form to the Lagrange-Rayleigh ID strain dynamics of Bales and Gooding (BG), with 'strain accel- 
erations' proportional to a Laplacian acting on a sum of the free energy strain derivative and frictional strain 
force. The tensorial St. Venant's elastic compatibility constraints that forbid defects, are used to determine the n 
non-order-parameter strains in terms of the OP strains, generating anisotropic and long-range OP contributions 
to the free energy, friction and noise. The same OP equations are obtained by either varying the displacement 
vector components, or by varying the N strains subject to the N c compatibility constraints. A Fokker-Planck 
equation, based on the BG dynamics with noise terms, is set up. The BG dynamics corresponds to a set of 
nonidentical nonlinear (strain) oscillators labeled by wavevector k, with competing short- and long-range cou- 
plings. The oscillators have different 'strain-mass' densities p(k) ~ l/k 2 and dampings ~ 1/ p(k) ~ k 2 , so 
the lighter large-A: oscillators equilibrate first, corresponding to earlier formation of smaller-scale oriented tex- 
tures. This produces a sequential-scale scenario for post-quench nucleation, elastic patterning, and hierarchical 
growth. Neglecting inertial effects yields a late-time dynamics for identifying extremal free energy states, that 
is of the time-dependent Ginzburg-Landau form, with nonlocal, anisotropic Onsager coefficients, that become 
constants for special parameter values. We consider in detail the two-dimensional (2D) unit-cell transitions from 
a triangular to a centered rectangular lattice (N op = 2,n = 1,N C — 1); and from a square to a rectangular 
lattice (Nop = 1, n = 2, N = 1) for which the OP compatibility kernel is retarded in time, or frequency- 
dependent in Fourier space (in fact, acoustically resonant in ui/k). We present structural dynamics for all other 
2D symmetry-allowed ferroelastic transitions: the procedure is also applicable to the 3D case. Simulations 
of the BG evolution equations confirm the inherent richness of the static and dynamic texturings, including 
strain oscillations, domain-wall propagation at near sound speeds, grain-boundary motion, and nonlocal 'elastic 
photocopying' of imposed local stress patterns. 



PACS numbers: 64.70.Kb, 64.60.Cn, ll.10.Lm, 81.30.Kf 

I. INTRODUCTION 

Structural phase transitions in solids have attracted a great 
deal of interest over a century, both for their conceptual impor- 
tance as symmetry-changing phase transitions, and for their 
role in inducing technologically useful materials properties. 
Both the diffusion-controlled replacive and the diffusionless 
displacive transformations have been studied, although the 
former have received more attention because their reaction ki- 
netics is more conducive to control. 

We consider here the class of materials known as ferroe- 
lastic martensites. Ferroelasticity is defined by the existence 
of two or more stable orientation states of a crystal that cor- 
respond to different arrangements ottbe atoms, but are struc- 
turally identical or enantiomorphousutl. In addition, these ori- 
entation states are degenerate in energy in the absence of me- 
chanical stress. The term 'martensitic' refers to a diffusionless 
first order phase transition which can be described in terms of 
one (or several successive) shear deformation(s) from a par- 
ent to a product phaseE The morphology and kinetics of the 
transition are dominated by the strain energy. The transition 
results in a characteristic lamellar or twinned microstrcuture. 

Salient features of ferroelastic crystals include mechani- 
cal hysteresis and mechanically (reversibly) switchable do- 
main patterns. Usually ferroelasticity occurs as a result of a 



phase transition from a non-ferroelastic high-symmetry 'par- 
ent' phase and is associated with the softening of an elastic 
modulus with decreasing temperature or increasing pressure 
in the parent phase. Since the ferroelastic transition is nor- 
mally weakly first order, or second order, it can be described 
to a good approximation by the Landau theory with sponta- 
neous strain or deviation of a given ferroelastic orientation 
state from the parent phase as the order parameter. The strain 
can be coupled to other fields such as electric polarization and 
magnetic moment and thus the crystal can have more than one 
transition. Depending on whether the spontaneous strain is the 
primary or a secondary order parameter at a given transition, 
the lower symmetry phase is called a proper or an improper 
ferroelastic, respectively. While martensites are proper fer- 
roelastics, examples of improper ferroelastics include ferro- 
electrics and magnetoelastics. 

There is a further subset of ferroelastic martensites (either 
non-elemental .metals or alloy systems) that exhibit the shape 
memory effecto. These materials are characterized by highly 
mobile twin boundaries and (often) show precursor structures 
above the transition. Furthermore, these materials have small 
Bain strain, elastic shear modulus softening, and a weakly to 
moderately first order transition. Some examples include InTl, 
FePd, NiTi and AuCd. 

Dynamics plays a central role in proper ferroelastic 
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transitions! 
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As noted, these materi- 
als undergo diffusionless, displacive transitions, with strain 
(components) as the order parameter, and develop complex 
microstructures in their dynamical evolution, finally form- 
ing spatially varying, multiscale 'textures' or strain pat- 
terns. When quenched, some martensitic materials de- 
velop interfaces moving at near-sound speeds. Textured im- 
proper ferroelastics include materials of .technological impor- 
tance such as superconducting cupratesEl and colossal mag- 
netoresistance (CMR) manganitesLj. Many dynamical mod- 
els have been invoked ^foJ]|jQW|-aM^ ferroe- 
lastic pattern formationfflQ&BlllUELjoOt^ such as nu- 
cleated twiii-fcpnt propagation; width-length scaling of twin 
dimensionsE2E!l; tweedE3tJc3; stress .effects; elastic domain 
misfits; and acoustic noise generationEl n 
In a one-dimensional (ID) model, Bales and Goodinga 
considered a displacement u and a sixth order free energy 
F(e), nonlinear in the strain, that in ID is simply a deriva- 
tive, e = Jiu/dx. With Lagrangian dynamics, a Rayleigh 
dissipationEj and variation in u, a single strain evolution equa- 
tion was obtained in ID: 
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where A' and po are the scaled friction (coefficient) and mass 
density, respectively. In a low-frequency/large-wavevector 
regime p^oj « A'k 2 , where the inertial term is small, a 
simple timgpdependent Ginzburg-Landau (TDGL) equation 
is obtainedoQ in II?: 
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For the k ^ modes in ID the strain is a one-component vec- 
tor that in Fourier space is simply proportional to the scalar 
displacement; and there is only one type of lattice symmetry. 
In higher dimensions D > 1, the strain (displacement) is a 
tensor (vector), and there are many possible discrete lattice 
symmetries. The central question is: What is the general form 
of the underdamped evolution equations for (order-parameter) 
strain-tensor components, for ferroelastic transitions of differ- 
ent symmetries? 

A wide variety of dynamical models have been used to date. 
These include a 2D or 3D TDGL dynamics in morphologi- 
cal variables for structural variants, with a long-range poten- 
tial between squarep-pf these variablesia, motivated by ideas 
of elastic inclusions^]. Other work in 2D and 3D has used a 
TDGL equation for the OP strains only, with a longrfaage 
potential emerging between the OP strains themselvesM-M 
(not thejj- squares), by optimizing non-OP strains. Some 
authorsu assumed the validity of a 2D Bales-Gooding (BG) 
form from a Lagrangian without non-OP strains, and then 
phenomenologically added a long-range potential between 
squares of the OP strains—A TDGL equation in the displace- 
ments has also been usedtd Yet other models considered dis- 
placement accelerations equated to displacement-gradients of 
the free energy as forces; otstrain TDGL equations coupled 
to vacancy-field dynamics^. Finally the Lagrange-Rayleigh 
procedureEEcJ has recently been appliedtJL] in 2D, yielding 



an underdamped dynamics for the displacement, that is trun- 
cated to overdamped equations, that are seemingly different 
from the TDGL form. 

While these and other models yield valuable insights 
into ferroelastic texturing (i.e., single-crystal microstructure), 
there is clearly a need, through explicit derivation, to obtain 
an underdamped, symmetry-specific OP strain dynamics for 
D > 1; to find the precise form of long-range potentials (if 
any) that emerge; and to determine the regime of validity (if it 
exists) of some form of TDGL equations. 

In this paper, we use the Lagrange-Rayleigh variational 
procedure to derive a ferroelastic strain dynamics (including 
noise terms). A central j^ejA-Blayed by the N c St. Venant 
compatibility conditionsEfEaoEy for the N = N op + n sym- 
metry adapted strains, that enforce the absence of defects 
(and lattice integrity) at each instant, and allow the n non- 
OP strains to be expressed in terms of N op order parame- 
ter strains. We show that: (i) an underdamped set of N op 
equations can be obtained for the OP strains alone, that is of 
a generalized BG form, with naturally emerging anisotropic 
long-range (ALR) contributions to OP potentials, friction and 
noise. (For the N op — 1 case these are in general also ex- 
plicitly retarded in time.) (ii) The same OP equations can 
be obtained, either by varying the displacement, or by vary- 
ing the strains subject to the compatibility constraint through 
dynamic Lagrange multipliers, (iii) Dropping strain inertial 
terms yields strain TDGL equations, with nonlocal Onsager 
coefficients, that reduce to constants for special friction val- 
ues, resulting in a local TDGL dynamics. These act as a 
late-time dynamics for the damping envelope of textural os- 
cillations. 

We explicitly demonstrate (i), (ii), (iii) above for the 2D 
triangular to centered rectangular or TR lattice transition 
(N op = 2, n = 1, JV C = 1) and for the square to rectangu- 
lar or SR lattice transition (N op = 1, n = 2, N c = 1), as well 
as present dynamics for all other allowed 2D symmetries. The 
procedure can be generalized to 3D, e.g. the cubic to tetrag- 
onal (N op — 2, n — 4, N c — 6) transition where (static) 
compatibility potentials in a TDGL dynamics produce rich 
texturestil Our central result is a generalized BG dynamics 
written in the OP strains {si}, £ = 1,2, N op only: 
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(1.3) 



where q is a symmetry-specific constant, and p n is a scaled 
mass density. F c ({ei\), R c ({ei}) are the compatibility- 
induced symmetry-specific contributions that emerge natu- 
rally from the non-OP free energy as additions to the OP 
free energy F and OP Rayleigh dissipation R, while gl is 
the corresponding noise term that adds to the OP noise gi. 
In equation (1.3) and subsequently, we use the symbol A 
to denote dimensionless discrete derivatives on a reference 
latticed. The generalized BG equations can be written as 
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Langevin equation for £g(k, t) and V£ = £e(k, t), yielding sta- 
tistically equivalent Fokker-Planck equations for the probabil- 
ity P({si, v t ], t). 

In the strain-variation derivation of the BG dynamics 
above, we introduce the concept of a strain mass-density ten- 
sor whose components (in Fourier space) behave as p aa i (k) ~ 
Po/k 2 , that is responsible in coordinate space, for the Lapla- 
cian on the right-hand side of (1.3). This generalization of 
the ID case expresses the physical idea that long-wavelength 
strains are extended lattice deformations and hence have 
greater inertia. We present a physically illuminating analogy 
of (generalized) BG dynamics as an array of coupled nonlin- 
ear fc-space oscillators that have an intrinisically hierarchical 
equilibration, with large k oscillators damping out first. 

The N op strain order-parameter equations with derived 
anisotropic long-range terms are equivalent to the D displace- 
ment equations that do not explicitly have such terms. The ad- 
vantage of the OP strain approach is that it displays and uses 
such anisotropic long-range correlations, that are valuable in 
understanding simulated textures, as demonstrated below. 

More generally, in the displacement (u) picture, the strains 
are derived quantities and the compatibility condition is an in- 
cidental identity, expressing the single-valuedness of u. This 
is analogous to describing magnetic problems in terms of the 
vector potential A, with B just a label for A x A, and with 
A • (A x A) = just expressing a vector identity. By con- 
trast, in the strain-only picture, the strain tensor components 
E^ v are the physical variables, and the compatibility condi- 
tions A x ^A x E^j =0 (with T denoting transpose) are 

treated as independent field equations expressing the physical 
constraint of no defects. This is analagous to working with 
the magnetic induction B, where the Maxwell's field equation 
A • B = expresses the absence of magnetic monopoles. The 
compatibility equation for strain-has been used for a consistent 
description of forces in liquidsEj; here we use it to develop a 
consistent OP dynamics for ferroelastics. 

The basic idea is quite simple. Compatibility implies non- 
OP strains {e^} are proportional in Fourier space to OP 
strains (e^): 

e i = Y^ l Sue l . (1.4) 

i 

Hence, the non-OP free energy / (and Rayleigh dissipation), 
harmonic in the non-OP strains (and time derivatives), can be 
written in terms of OP strains. With elastic constants a^, 

/ = \ E^N 2 = \Ys Uuietee* = F c ({ £i }). (1.5a) 

i.k i,i'k 

This defines the compatibility kernels 

Uw = y^^jSuSw* (1.56) 

i 

for F c (and similarly for R c ) of (1.3) in the desired OP-only 
dynamics. Thus the problem reduces to finding the propor- 



tionality constants Su, for each symmetry -based phase transi- 
tion. 

The plan of the paper, with self-contained Sections, is as 
follows. In Sec. II the OP dynamics for the TR case is 
derived by u variation. In Sec. Ill we demonstrate that the 
same TR dynamics is obtained by strain variation, with en- 
forced compatibility. Results for the SR case are stated and 
analyzed, with derivation details in Appendix A. Numerical 
simulations of some interesting BG dynamic^volutions are 
presented (using standardized scaled energiesEI and dissipa- 
tions, as in Appendix B). Noise contributions are derived, 
and a Fokker-Planck formalism is set up. Section IV dis- 
cusses the equivalent inhomogeneous oscillator description of 
the BG dynamics. Section V deals with its TDGL trunca- 
tion, also .derived in Appendix C from truncated displacement 
dynamical§0. Section VI presents the compatibility kernels 
for other 2D symmetries. Finally, Sec. VII contains a sum- 
mary and discussion. 

II. ORDER PARAMETER STRAIN DYNAMICS BY 
DISPLACEMENT VARIATION 

Consider a Lagrangian density L(a,a) = J dtJ2?(T ~ 
V) that depends on a variable a(f, t) through a kinetic term 
T = T{a) = J2rlP°® 2 ' and ajjotential term V = V(a). 
Then, with a Rayleigh dissipatior£3 R tot = j dr^r/a 2 (where 
T] is the friction coefficient), we have, by variation in a, the 
Lagrange-Rayleigh equation: 

d dL dL _ dR tot 

dt da da da 

In this Section we work in the displacement picture, and 
consider variations in displacement a(r,t) — ► u M (r, i) that 
will generate \i = 1, ,.,D equations, for it p in D dimensions. 
The potential or Gibbs free energy V depends on displace- 
ment derivatives, that are = N op + n symmetry-adapted 
linear combinations of the symmetric Cauchy strain tensor 

E„ v = i (A M u„ + A„ Ull ) . (2.2) 

We neglect 'geometric' nonlinearity, as justified by the scal- 
ing of Appendix B, where such corrections are higher order 
in a typical (small) strain value. We notationally distinguish 
between OP strains {e^} and non-OP strains {e^}. The po- 
tential V = F({ee}) + /({ej}) is anharmonic through F 
in the N op order parameter strains, and harmonic through 
/ in the n non-OP strains, while the Rayleigh dissipation 
R tct ({£e}, {ii}) is harmonic in both the strain rates: 

r,i r,£ r,i 

Here as in Appendix B, {Ag\, {a^} and {A' e }, {a'A- are, 
respectively, OP and non-OP second order elastic and fric- 
tion coefficients, and the sum is over sites {r} of a reference 
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lattice, while t is a scaled time. 

TR dynamics from displacement variation 

Consider the triangular to centred rectangular lattice (TR) 
transition, for which N op = 2 and n = 1. Figure 1 shows 
the TR, the square to rectangular or SR, and other lattice 
transitions. (While it is trueHtZl that these correspond to 2D 
projections/analogs of hexagonal to orthorhombic, and tetrag- 
onal to orthorhombic lattice transitions, respectively, we will 
reserve this 3D terminology for full 3D analyses, elsewhere.) 
The symmetry-adapted non-OP compressional strain ei = 
i (A x u x + AyUy), whereas the OP are the 'deviatoric' £2 = 
i (A x u x — AyUy) and shear strain £3 = i [A x u y + A y u x ). 
The deviatoric strain can be regarded as a diagonal shear. 

Then (2.3) for the non-OP compressional energy is f = 
£4«i e i 2 > R tot = |EfKe? + (A' 2 i 2 + A' 3 e\)]. The an- 
harmonic and fourth-order (triple-well) free energy F for £ 2 , 
£3, is given below in (3.27a), although this explicit form is not 
needed in the derivation. 

Defining OP free energy derivatives (i.e., stresses) ^2,3 = 
Se2s{r t) ' Lagrange-Rayleigh variation with respect to dis- 
placements u(r, t) gives for the dynamics of (2.1) 



1 dF 
PoS2 = -^{ai(A x 2 - A y 2 ) ei + A 2 ^-] 

+ i[al(A* 2 - Ay 2 )^ + A' 2 A 2 e 2 }; (2.6b) 

P0S3 = -[2a 1 A x Aye 1 +A 2 —] 
4 oe 3 

+ ^[2a' 1 A x A y e 1 + A' 3 A 2 e 3 \. (2.6c) 

By taking appropriate derivatives of (2.6), it is easy to see 
that the compatibility condition (2.5) is satisfied as an identity. 
This linear equation, (2.5), can then be used instead of the 
nonlinear (2.6a), to eliminate the non-OP strain ei(k,oj) in 
terms of the OP strains e 2i 3(k, lo), assuming periodic bound- 
ary conditions. 

Defining ai u = a\ — ioja'x, transforming to Fourier space, 
and with P2 3 now defined as P2 3 = — — , the compati- 

bility constraint (2.5) 



Pou x = -[aiA^ei + A X F 2 + A y F 3 ] 



-KA^ei + A' 2 A x i 2 + A' 3 A y i 3 ]; (2.4a) 



Potiy = -[aiA. y ei - A y F 2 + A X F 3 ] 



+ -[a[A y e 1 - A' 2 A y e 2 + A' 3 A x e 3 }. (2.46) 

These displacement equations have been obtained 
previously^, but were then truncated by dropping the 
displacement acceleration (second time derivative) to yield 
reduced equations that are analyzed further in Appendix 
C. Instead, we pursue here the underdamped OP-strain 
equations and find they have a generalized BG form. 

The strains obey the compatibility constraints, that in the 
displacement picture ensure that u is single-valued (i.e., cross- 
derivatives commute). The equation for the 2D case is, at 
every instant, 



A 2 e: 



A 2 ) £ 2 - 2A.A.E3 = 0. 



(2.5) 



Taking spatial derivatives of (2.4) we obtain the full under- 
damped equations for the strains: 

1 -* dF dF 

Poe'i = 7 [aiA 2 ei + (A x 2 - A y 2 ) — + 2A x A y —\ 
4 oe 2 063 



-KA 2 ei + A' 2 (A X 2 - A 2 )e 2 + 2A' 3 A x A v i 3 }; (2.6a) 



Q iei (k,u) + Q^ik^u) + Q 3 e 3 (k,u) = (2.7) 

with Qi = — fc 2 , Q 2 = k x 2 — k y 2 , Q 3 = 2k x k y , is used to 
eliminate e\(k-u), giving Nop = 2 order parameter strain- 
only equationsEH: 



PqUJ £ 2 = — 



PqUJ £3 = — 



Q\ Q2Q3 
Qf 2 Qi 



£3 ) + F 2 - iLoA' 2 e 2 
(2.8 



Q2Q3 . Q 2 



Qi Qi 



£2 + ^§£3 + F 3 - iuA' 3 e 3 



(2.86) 

These can be succinctly written as an OP dynamics for 
e t (r, t) with (£ = 2,3): 



p et = -A 



1*2 f8(F + F c ) , S(R + R C ) 



See 



See 



(2.9) 



where F c and R c are compatibility-induced contributions 
from the non-OP free energy / and dissipation R written in 
terms of the OP strains: 



F c = J2 U U' &)ei{k, t)e* t > (k, *), (2.10a) 

k,e,e> 



RC = 2 a 'i E Ju>$)£t$,t)<?i'$,t), (2.106) 

where the orientation-dependent kernel0 U[ e ,(k) ~ Vu'{k) is 
defined implicitly above, and given explicitly in (3.29) of Sec. 
Ill, following a strain-based derivation of the same dynamics. 
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The system of N op = 2 underdamped equations (2.9) 
derived for the strains is clearly of a Bales-Gooding form 
[compare with (1.1)], but now generalized in three ways: 
by the derived replacement -^i — ► A 2 ; by the appearance 
of a compatibility-induced anisotropic long-range (ALR) in- 
teraction, between ei-ey (and not e 2 -e 2 ,); and of a similiar 
compatibility-induced ALR dissipation, between ig-if,' strain 
rates. 

An interesting consequence of the BG structure with pe- 
riodic boundary conditions is that the "fc = 0" (or more pre- 
cisely, k — > 0) OP strain obeys poi'e(k = 0, t) = 0, so macro- 
scopic strain momentum poSt(k = 0, t) — is conserved^], 
unaffected by internal forces and dissipations. The solution 
is eg(k = 0,t) = et(k = 0,0)t + e e {k = 0,0). In 
the special case of an austenite phase with initial conditions 

ii{k = 0,0) = 0, and e e (k = 0,0) - Er £ *( r >°) = °> 
strains of both signs will develop on cooling below transition, 
as a consequence of the dynamics. The OP strain is in this 
case like a 'charge' that is generated in sign-balancing pairs, 
and the notion of elastic 'screening' helps in understanding 
simulations presented later. 

The square-rectangle (SR) dynamics, driven by a devia- 
toric strain OP, can similiarly be shown by displacement vari- 
ation to be also of the generalized BG form, as given in the 
first part of Appendix A. However, we now proceed to derive 
the strain dynamics through strain variation. 



III. OP STRAIN DYNAMICS BY STRAIN VARIATION 
WITH COMPATIBILITY CONSTRAINTS 

In discussions of ferroelastics, it is common to assert that 
although the free energy is in terms of the strains, the true ba- 
sic variables for such systems are displacements, since strains 
are just displacement derivatives. Thus Monte Carlo simu- 
lations, numerical solutions of dynamic equations and static 
analyses of textures, even when expressed in terms of strains, 
are all finally performed in terms of displacements. Following 
the electromagnetic analogy mentioned in the Introduction, an 
alternative treatment is in terms of strains as the basic vari- 
ables. The free energy for a first order transition, say, Fo ~ e 6 , 
is then regarded as zeroth order in derivatives, whereas in the 
displacement picture it is a sixth power of derivatives. In 
this Section, we derive ferroelastic dynamics for the TR and 
SR cases, using strains as the variational quantities. Results 
for external stress and noise are also_siated. We use peri- 



odic boundary conditions throughouO^J and transform be- 
tween coordinate r and wavevector k descriptions, as conve- 
nient. The Lagrange-Rayleigh dynamics equations for a gen- 
eral variable a is 



d dL 
dt da 



8L 
da 



dR l 



da 



and we consider here the strains as the variables: 



a^{E^(f,t)}. 



(3-1) 



(3.2) 
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FIG. 1: Symmetry-allowed transitions in 2D for the four crystal sys- 
tems. The dark lines are guides to the eye, for deformations. There 
is a one-component strain order parameter for: (a) the square to rect- 
angle or SR case, driven by £2; (b) the rectangle to oblique or RO 
case driven by £2; and (c) the square to centered rectangle or SC 
case, driven by £3. A two-component OP, or two one-component 
OP's, lead(s) to (d) the triangular to centered rectangle or TR case, 
driven by £2, £3; (e) the triangle to oblique or TO case, driven by 
£2, £3; and (f) the square to oblique or SO case, driven by £2 and 
£3, independently. [The oblique to oblique ferroelastic (PI to PI) 
transition that involves merely a loss of inversion symmetry is not 
consideredu.] 



In total, if we consider.-all space groups in 2D, there are 
23 ferroelastic transitionstJ in 2D. Figure 1 shows the six 
symmetry-allowed transitions in 2D for the four crystal sys- 
tems. However, of these only six correspond to mono-atomic 
basis. We will present the initial part of the argument in gen- 
eral form, before focusing on the triangular to (centered) rect- 
angular lattice or TR case; and the square to rectangle or SR 
case, with other 2D transitions considered in Sec. VI. 

The Lagrangian contains the potential V and kinetic energy 
T that depend on symmetry-adapted strains or linear combina- 
tions of the strain tensor {E^}, that are irreducible represen- 
tations of the unit-cell symmetry group. The s = 1, 2, . . . , N 
symmetry-adapted strains are written as {e s }, with N = 3 in 
2D (and 6 in 3D). In 2D the compressional (ei), deviatoric 
(e2) and shear (es) strains are defined by 



£1 
ci 



1 



■(E XX +Eyy), 



£2 

C2 



1 



c 3 2 

where c\, C2, C3 are symmetry-specific constants. 



(3.3) 
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The strain tensor (neglecting geometric nonlinearity as jus- 
tified in Appendix B) obeys the St. Venant compatibility con- 
dition, that is here a field equation forbidding defects such as 



dislocations and vacanciestxi 



at every instant: 

T 



A x A x E(r*, t) = 0, 



k x E(k, t) x k = 0, 



(3.4a) 



(3.46) 



with no source term on the right hand side. The p 
1,2, . . . , N c compatibility conditions are: 



(3.5) 



where Qji are second-order derivative operators, from (3.4a). 
In 2D there is only one compatibility equation N c = 1, and 
the operators in Fourier space are Qi (k) = —k 2 /ci,Q 2 (k) = 
(k x 2 - k y 2 )/c 2 , Q 3 {k) = 2k x k y /c 3 . Thus (3.5) for e s (k,t) 
is, from (3.4b), 

Qiei + Q 2 e 2 + Q 3 e 3 = 0. (3.6) 

The anisotropic compatibility factors Qi j2) 3(fc) encode the 
discrete symmetries of the compression, deviatoric and shear 
strains. The symmetry constants for the TR case are c\ = 
Co = C3 = 1, whereas for the SR case they are c\ — c 2 = 
V2, C3 = 1, The compatibility constraint will be invoked re- 
peatedly in the derivations below. The physical meaning of 
the constraint is that order parameter strains should not tear, 
or cause defects in the lattice, i.e. lattice integrity is main- 
tained. Suppose, in a sea of square unit cells, one cell was 
made rectangular (local deviatoric strain). It is clear that in 
order to maintain lattice integrity, the neighboring cells must 
also deform, inducing all three strains in an inter-related way, 
with a similar (but smaller) deformation of the larger number 
of further neighbors. The requirement of a smooth, compat- 
ible fitting together of neighboring unit cells will cause the 
disturbance to propagate outwards, and in an anisotropic way 
(due to discrete crystal symmetry and elastic constants): the 
local condition has global consequences. 

The Lagrangian L also contains the compatibility con- 
straints through dynamic Lagrange multipliers {A^(f, t)}: 



L 



dt 



(3.7) 



We now: (a) obtain the kinetic energy T({e^}) in terms of 
the time derivatives of the symmetry-adapted strains; and (b) 
use this to derive the TR and SR dynamics by strain variation, 
incorporating the compatibility constraint. 



A. Kinetic energy in terms of strain rates: 

Since Newtonian dynamics is for point particles, the kinetic 
energy in terms of displacements is: 

T = ^° E *) = kf* E M*> *)i 2 ( 3 - 8 ) 



in coordinate and wave-vector spaces, where po is a dimen- 
sionless mass-density that is a ratio of typical kinetic and 
elastic energy densities (Appendix B). In the strain picture, 
the displacement can be derived fcpm the strains through the 
Kirchoff-Cesaro- Volterra relationEI 



*)= f [E(f) + {(I- r) x V f x E(f)}].dT, (3.9a) 

JC(r ,r) 



where the line integral is along any contour C(f , t) to f from 
f , that is a fixed point of the deformation. Taking derivatives 
with respect to r, the symmetric combination of derivatives is 



h^u v + A v u„) = E lxv . (3.96) 

In 2D, using (3.3) this gives 

AjMj = 1 ; A v u u = . (3.10a) 

Cl c 2 Ci c 2 

Taking time derivatives and transforming to Fourier space, 



Ux(k, t) = 



_|_ e2_ ei_ e% 



Uy(k,t) 



Ik,, 



(3.106) 



Inserting (3.10b) into (3.8) yields the kinetic energy, that is 
nonlocal in terms of the strain rates: 



T = E \pss'{k)e*{k,t)e s ,{k,t) 



E \psAr-?)eW,t)eA?,t), (3.11) 



where we have introduced an anisotropic 'strain mass-density 
tensor' whose components turn out to be related to ratios of 
the compatibility factors of (3.6): 



p ss >(k) = p{k) 



1 fe 2 

1 (*;-*; 

cic 2 



cic 2 fep| 



C 2 ^x^y 



P (k) 



(Qi\2 ( QiQi \ 
( Q2Q1 \ c 2 i(Ql\2 



(3.12a) 



This strain mass density tensor is a kinematic time- 
independent quantity true for all 2D symmetries and has a 
determinant (2 / ' ciC2k x k y ) 2 . Here 



p{k) 



4p 



(3.126) 



r.ji 



and therefore long-wavelength strains over many lattice spac- 
ings are effectively more 'massive', as is physically reason- 
able. In 2D coordinate space, p ss i (R) ~ (Ap )ln(\R\). It is 
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this inverse Laplacian dependence that gives rise to the Bales- 
Gooding structure (1.3) of the underdamped dynamics. The 
strain kinetic energy can be expressed only in terms of the 
compressional and deviatoric strain rates, i.e., the 'shear com- 
ponents' of the strain mass tensor are zero, p S 3 = ps s = 0. 
This is because from (3.3) and (3.9b), the shear strain rate 

= \[ik x u y {k) + ik y u x (k% (3.13) 
C3 ^ 

is not independent, but is related to the other strain rates by 
a consistency condition through (3.10b), that turns out to be 
precisely the compatibility constraint. 



B. Dynamics by strain variations: 

It is convenient to henceforth notationally distinguish be- 
tween n non-order parameter strains {e^} and N op order pa- 
rameter strains {eg}. The compatibility conditions (3.5) be- 
come 

C (p) (r, i) = £ Q { e p) ee(r, t) + E Q^e^r, t) = 0, 

£ i 

(3.14a) 

and therefore in Fourier space 

C<*)$, i) = E Q^\k)et(k, t) + J2 Q\ p) (k) ei (k, t) = 0. 

l i 

(3.146) 

The Gibbs free energy V({ei}, {e^}) depends harmonically 
on the non-OP strains through / and anharmonically on OP 
strains through F, whereas the Rayleigh dissipation R tot to 
the lowest order depends harmonically on all strain rates. 
Thus, 

V = f({e i })+F({e i }); / = ^ 5>e?; 

R tot = lj2^i 2 + lYl A ^l ( 3 - 15 ) 

where o», a[ are the non-OP elastic and friction constants. 

The kinetic energy in terms of strain rates {e'j}, {ie}, from 
(3.11), is 

T = E \pw(k)e* t {k,t)e t (k,t) 

U,t,f 

+ _E \piAk)e*{k,t)e v {k,t)}. (3.16) 

Using (3.14b), (3.15) and (3.16), the Lagrange-Rayleigh dy- 
namics for OP and non-OP strains are given by directly vary- 
ing the strains in (3.1) and (3.7): 

dF 



E Pu >h> + E = - E Qt )A<p) ~ 

p 

(3.17a) 



E P«' e 'i' + aiei + + E = ~ E PU'£&' > 

i' p I' 

(3.176) 



E 



E^W 



(3.17c) 



We have written the equations in a general form for future 
3D BG .-generalizations such as the cubic to tetragonal 
transitionllil. There are (n + N c ) linear equations (3.17b,c) for 
(n + N c ) variables {e^}, {A' p '} that can be written in matrix 
form, and inverted to yield the dynamics for the OP strains 
{ei}. We do not pursue this general treatment here, but now 
specialize to the TR case and give the result for SR at the 
end of the section. Other symmetries are considered in Sec. 
VI. 

TR underdamped dynamics: 

For the TR transition, n = 1, N op = 2, N c = 1; the 
non-OP strain is compressional e\\ the OP strains are 'de- 
viatoric' and shear (e 2 , £3 respectively); whereas the symme- 
try constants of (3.6) are c\ — c 2 = C3 = 1, The strain 
mass-density components are pn = p{k)(Qi/Qj,) 2 = P22, 
pia = p(k)(QiQ 2 /Q3 2 ) = P2i- We have N op + 
2 + 1 + 1 = 4 equations like (3.17): 



Nr. 



Pn£l + Pi\&\ 



dF 

ds 



dF 



- - O2A - A' 2 i 2 ; 



3 A- A'e 3 ; 



(3.18a) 



(3.186) 



Pne'i + pi2£2 = —aid - a[ei - QiA; (3.18c) 



Qiei 



- Q,3£3- 



(3.18d) 



There are N. 



op 



2 nonlinear equations for e 2; 3 and 
n+N c = 2 linear equations for e±, A. We can use the compat- 
ibility equation (3.18d) for the strain mass density {p ss i(k)} 
components in (3.12) to write the expression on the LHS of 
(3.18c) as 



.. . .. 4p .. 

pliei + pi2£2 = 77- -£3- 
W3 



(3.19) 



The Lagrange multiplier of (3.18c) is determined by the 
OP strains as 



A(fc,w) = ^ 



ai^£2 + a luJ Q 3 + 



Q3 



£3 



(3.20) 

where a\ u = a\ — iu>[. Substituting (3.20) into (3.18) with 
the identity 



O2 4^ 

P22 - P2I7— " — 7 -jr, 

Qi k l 



(3.21) 
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yields the TR equations in terms of the OP strains alonel 



El 



Pquj e 2 = — 



p uj 2 e a = 



Se* 2 
SF_ 



Qt 

Q 2 i 



Q2Q3 

— —e 3 J - iujA 2 e 2 



(3.22a) 



Q2Q3 , Q 



£2 + t4 £ 3 - iuA'a£ 3 



Qi ' Qi 

(3.226) 

Since ratios of compatibility factors recur, it is useful to 
define 



Qt,t{k) 



Qt(k) 



Qr(fc) 

The non-OP strain in (3.18d) is a derived quantity 

ei(fc,u>) = - Q,2,i(k)e 2 (k,uj) + Q 3 ^{k)e 3 {k,u) 



(3.23) 



The compatibility induced OP-OP interaction is the 
non-OP free energy written in terms of the OP strains, 

/({ei({e£})})^F c ({££}) with (£ = 2,3): 



n{f,t)Ui l ,{?-i*)e l ,{i>,t) 



UiY, U ( c j,Ck)s e (k,t)e* e ,(k,t). (3.276) 



The Rayleigh dissipation function for the OP is 



1 



r,(. 



(3.28a) 



and the compatibility-induced contribution R c is the non-OP 
dissipation written in terms of the OP strain rates 



S12E2 + S 



13=3) 



(3.24) 



defining for the TP case, the constants Su(k) mentioned in 
(1.4), that depend only pn the wavevector direction k, i.e., no 
non-trivial length scaleEil, and not on its magnitude |fc| (i.e., 
no length scale). 

In a compact form with i = 2,3, 



2 k 2 / 5{F + F c ) 5{R + R c ) 

Pqlu eg = — - 1 = 

4 V Se*(k,uj) 5iUk,u) 



(3.25) 



This is written as an OP strain-only dynamics for ee(r, t) 



.. , . 1., M(F + F C ) S(R + R C )\ ,„„„, 
Po^r.t =-A 2 i / + \, ' , 3.26 

that can-, be written as a strain-momentum continuity 
equation!!! The OP free energy from Appendix B is 



+{(e 2 2 + £ 3 2 ) - 2(4 - 3e 2 el) + [e\ + el) 2 }]. (3.27a) 

This free energy is invariant under the 6mm point group oper- 
ations. Specifically, under three fold rotations, E x *x f = -%^ + 



3Ey 



~s/3 p . 771 
2 ^xyi ^y'y' 



1 TP 

4 -^ica; 



1 171 _ \/3 171 . 771 



— ^{E xx — E yy ) — |P K y an( J the non-OP and OP strains 

transform as e' x = e\, e 2 = —\^2 — "^£3> £3 = -^£2 — ^£3- 
The (anharmonic) OP free energy has not been explicitly used 
in the derivation of the dynamics, whose structure depends 
only on the number and nature of the (harmonic) non-OP 
strains. 



RC = T E k{r,t)r]lAr-n^,t) 



r, ?'£,£' 



} E ril e {k)u 2 h{k,t)e},{k,t). 



Thus, explicitly 

8{R + R C ) ^ 7 . 

TTTTf - =E A W fc ^)' 
dee \k,u) J, 

where the effective OP friction is 

A e.e> = K^i.i' + a 'iV c w 



(3.286) 



(3.28c) 



(3.28d) 



Here the (frequency-independent) potential and friction ker- 
nels emerging from the dynamics are the same, 

rf u ,[k)= U c w(k), (3.29a) 

where 

U c a' = SuSw = QisQi',i = Q1Q1 1 /Q\ 2 , (3.296) 
and (as from static-constrained minimization!!), 



= (Q2/0l) 2 ; 



^) = 4^ = (0 3 /0i) 2 ; 



A: 4 



(3.29c) 
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Both the friction and potential kernels depend on the waveji- 
tor direction k, are independent of \k\ for long wavelengthsM, 
and are 'anisotropic long-range' functions encoding the dis- 
crete symmetry of the lattice, that fall off in coordinate space 
with a dimensional power, r] c (R), U C (R) ~ l/R D . At k = 0, 
the potentials are undefined and we set U[ eU , = 0. (A 
Coulomb potential, by contrast, diverges for long wavelengths 
as - 1/k 2 and falls off as - 1/R D -' 2 .) 

Although (3.26), derived from strain variation will give the 
same results as (2.4), derived from displacement variation, 
they are conceptually distinct. In the displacement picture, 
(2.8), or any equation obtained from it, is solved with (u x , u y ) 
on a lattice, with both initial and boundary conditions applied 
to the basic variables u. Strains are defined as derivatives of 
the basic variables, and are derived quantities, e.g., F ~ ee* 
is a fourth power of derivatives. By contrast, (3.26) is an OP- 
strain-only dynamics, and is solved in the strain picture, with 
£2, £3 on a lattice; with initial and boundary conditions ap- 
plied to the basic variables {ei}\ and with the F ~ £( i term as 
zeroth order in derivatives. (Indeed, the highest order deriva- 
tives are strain-gradient squared or Ginzburg terms.) The non- 
OP strain is a derived quantity, obtained after solving for the 
OP and then using (3.24), the displacement can also, in prin- 
ciple, be derived using (3.9). 

In displacement-picture simulations, terms of the Landau 
free energy and Rayleigh dissipation are all anisotropic, 
being powers of the various displacement-derivative com- 
binations. In strain-picture simulations, the anisotropy is in 
the symmetry-specific compatibility kernels, with OP strains 
numerically treated as isotropic 'scalars'. The strain picture 
has advantages, as it works directly with the physical strain 
variables, and uses compatibility kernels evaluated once and 
for all, that encode the unit-cell symmetries and give insight 
into energetically favored textures. 

SR underdamped dynamics: 

We now turn to the square-rectangle or SR case, that shows 
a different structure, namely time-retarded OP potentials and 
friction. For the SR case , n = 2, N op = 1, N c = 1; the non- 
OP strains are compression and shear e±, e 3 ; the OP strain 
is deviatoric £2, the symmetry constants are c\ = c 2 = \/2, 
c 3 = 1. The compatibility factors are Q\{k) = — fc 2 /\/2, 



Q2 = {k% 
constraint (3.6) become 



k 2 v )/V2, 



2k x k y . Then the compatibility 



k 2 ex 



VSk x k y e 3 - (k 2 . - k 2 )e 2 



0. 



(3.30) 



The strain mass components are p\\ = p(k)(Qi/Q 3 ) 2 = /? 22 , 
P12 = p(k)QiQ 2 /Q3 2 = P2i- We have n + N op + N c = 2 + 
1 + 1=4 equations like (3.17), and with arguments similiar 
to the TR case we have, from the second part of Appendix A, 
the equations (A9) for £2, ei, S3: 



{puh + P12E2} = -oiei - a[ei - QiA, (3.31c) 



hei + Q,3e 3 = -Q 2 e 2 - 



(3.31d) 



Substituting for the compressional strain e\ yields coupled 
equations between e 3 and €2- As in (3.23) it is convenient to 
define a variable for the ubiquitous Q ratios, through Qe,e' = 
Ql/Qt 1 - Then, 



P0£2 



C2 2 k 2 



?3,i( fl i e 3 + a ie 3 ) 



(01QI i£ 2 + F 2 ) + KQ 2 ,i 2 + A' 2 )S2 }, (3.32a) 



c 3 2 fc 5 



[<92,lQ3,l( a l £ 2 + ai£ 2 ) 



+ (a 3 + axQg i)e 3 + (a' 3 + a' 1 Q 31 )e 3 



(3.325) 



Eliminating e 3 , these yield the final result (A15) in terms of 
the OP strain alone: 



p u £2 



c 2 2 k 2 \6F , b u {Q 2 /Q 3 f 



Set 



£2 — iuM 2 £ 2 



(3.33) 

As 6 W = [a 3uJ - (Apo^/clk 2 )}/^^, where a lu = a t - wa' it 
the kernel is now frequency dependent, and complex, as is the 
connection (A13) to non-OP strains: 



ei — S\2£2\ Su 



63 — S 3 2£2', S 3 2 



(QiQ 2 /Q 3 2 ) j 



(Q2/Q3) 



(3.34a) 



(3.346) 



where B u = 1 + b u {^) 2 . The compatibility condition is 
manifestly satisfied as an identity. 

It is convenient to write the second term on the right-hand 
side of (3.33) in terms of real and imaginary parts 

^^^^^(l) 2 )-^,^,^) 2 ), (3.35) 
with U c , r\ c given explicitly later. With this separation, 

6(F + F C ) 5{R + R C ) 



C 2 2 ,2 



p u> e 2 = —k 



fc5(/c,w) 5£2(k,u>) 



(3.36) 



p 22 e 2 + p 2 iei 



SF_ 

Set 



Q 2 A - A' 2 i 2 , (3.31a) where c 2 2 = 2. The OP dynamics for £ 2 {r, t) is then 



= -a 3 e 3 - a 3 e 3 - Q 3 A, 



(3.316) 



Po£ 2 (r,t) = -A 2 



d{F + F c ) d{R + R c ) 



de 2 {r,t) di 2 (f,t) 



, (3.37) 
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where the OP triple-well free energy is the same as in Ap- 
pendix B: 

F = £>- l)e 2 2 + £2 2 (e 2 2 l) 2 + ^(A £2 ) 2 . (3.38a) 



7 7 c (fc,c 2 , W 2 ) = T =1 -=-, (3.40c) 



(l + &rQl,3 2 ) 2 +^(0l,3) 



Under four-fold rotations, F T 



E xy , the non-OP and OP strains transform as 



xy 

e\ — > ei;e 2 — * —£2; e.3 — * — e3, leaving the free energy in- 
variant and similarly for other 4mm point group operations. 
The (anharmonic) OP free energy Fq is not explicitly used in 
the derivation, and the dynamics depends only on the number 
and type of the non-OP strains. 

We note that ui dependences carry an infinitesimal imagi- 
nary part to maintain causality, and thus (3.37) has contribu- 
tions only from earlier times: 



dF c 
de 2 



f dt'^U c {r-r' 1 t-t')e 2 {r' 1 t'). (3.386) 

J — OO 



The OP dissipation in addition to 



(3.39a) 



now also has a retarded compatibility contribution, 



dR c 
de 2 



dt' ^^(f-r^i- t'). (3.396) 



Both are ALR functions U c , i] c ~ l/R D as before but are 
also retarded in time. For negligible non-OP friction a[ = 
a' 3 = 0, the compatibility-induced friction vanishes, rf = 0. 
The non-OP compressional and shear strains ei^{r,t) are 
derived quantities obtained after solving for e 2 (r, t) and us- 
ing (3.34), and the constants Su are frequency-dependent and 
therefore retarded in time, like the compatibility potential. 

The real and imaginary parts of the kernel of (3.35) that 
give rise to the compatibility-induced OP potential U c and 
the friction coefficient if are explicitly 



a 3 . 



*7 c (fcW) = -0 



1 + SI & "IW 



oi ' [(l + 6 r 0i,3 2 ) 2 + 6 2 0i,3 4 ]' 

(3.40a) 

where the 'velocity' v = uu/k, and 



^03(5+^11 

" ai (l + (^)2) 



, bi 



» ai a 3 



|6 w | 2 = 6 r 2 + 6 l 2 = 
Similarly 



[i + (-fr) 2 ] 

(V + (^) 2 ) 
(i + (^) 2 ) ' 



(3.406) 



with 



4p 
c|a 3 



(3.40d) 



Note that both the static and dynamic U c are zero for di- 
agonal orientations, when Q 2 = 0. The zero-frequency 
limit U c Lk-D,0) = U c (k) is the static bulk .compatibility 
potentiaOcj used in earlier TDGL simulations!!!], and is 



f/ c o(fc) = 02,i 2 /[l + (ai/a 3 )03, 



(3.41a) 



This favors diagonal k x /k y = ±1 textures (when O2 = 0), 
and the derived non-OP strains ei, e 2 are expelled in a kind 
of 'elastic Meissner effect'l!!], as can be seen from (3.34). 
The bulk static compatibility potential is independent of \k\ 
for long wavelengths and does not set a domain-wall separa- 
tion length scale. Near the center point (T) of the Brillouin 
Zone (BZ) the U§(k) depends only on the wavevector direc- 
tion. It is-pnly in combination with the surface compatibility 
potential^ u surface (k) ~ 1/k thalpjjual-width 'true' twins 
(that satisfy a width-length scalingE3E!l) are obtained!!]. 

An exact Fourier transform to coordinate spaceEj confirms 
this preference for diagonal orientations: with cos0 = f.r', 
the coordinate space compatibility potential 



U c (r-r>) = 



G(0) 



(3.416) 



is found to have in the prefactor, a basic four-lobe structure 
from cos 4$, with higher harmonics. With a = — 



G(B) = 



[5 cos 40 (1 



sin 26 ' 



(17 + 4^) sin 20 + (13 - 4^ 

v l-\-a z ' ^ l-\-a z 



16v^T+^)(l 



l+a- 



COS 20) 2 

(3.41c) 

Figure 9 in Appendix C (TDGL dynamics), with a single-site 
initial condition, clearly shows multi-lobe strain textures from 
these higher harmonics. 

Returning to the dynamic case, we see from (3.40) that 
the repulsive kernel has a velocity-resonant structure through 
b r ~ g = 1 — (v/vq) 2 , strongest at a propagation speed 
V = vq = \J a^/Apo, or a time-dependent propagation length 
scale L p (i) = t/t p with t p — (ipo/as) 1 / 2 . These non- 
OP inertial effects, with anisotropic directional modulation or 
finite-velocity retardation ~ (u/k) 2 , compete with frictional 



delays 



with the peaks becoming singular for vanishing 



friction. This is somewhat like electrons interacting by ex- 
changing a photon, where the finiteness of the speed of light 
produces retarded Coulomb potentials!! The resonant struc- 
ture can be thought of as the OP strain textures 'exchanging 
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a non-0 P phonon mode', causing inertial time delays. L p (t) 
is the time-dependent size of an expanding anisotropic 're- 
gion of influence' within which changes in texturing at one 
point enforce compatible changes in texturing at other points 
of the lattice. This phonon mediation is also described by the 
equivalent instantaneous dynamics of (3.32) with deviatoric 
and shear strains obeying two coupled underdamped equa- 
tions, that are convenient for numerical simulationsH. 



nario for nucleation and growth. 





FIG. 2: Compatibility potential U c (k,uj 2 , (ui/k) 2 ) for the square- 
rectangle SR case versus scaled velocity, v = (ipo/as) 5 (uj/k), 
and direction parameter a — k x /k y for frequency/dissipation pa- 
rameter lj 2 = (ujd /a) 2 — 0.1. Thus U c is a repulsive and dynamic 
orientational potential, favoring diagonal textures, and moving inter- 
faces. 

The anisotropic potential kernel, U c (k,uj 2 , ) 2 ), is plot- 
ted in Fig. 2. This shows plots of the kernel versus scaled 

velocity, namely v = \J v and various BZ wave- 
vector directionsEl a = k x /k y for scaled frequencies Hi 2 = 
uj[(a! /a)] 2 = 0.1, where we take a[/ai = a' 3 /a3 = a' /a. U c 
determines the positive energy costs of the non-OP strains, 
that the system wants to eliminate. The smallest U c is at 
diagonal texturing |a| = 1, and zero velocity, v = 0. 
The most-non-optimal strains are the most strongly driven, 
with large U c . The striking features of Fig. 1 are clearly 
a preferred diagonal orientation a = ±1, and a strongest- 
repulsion textural velocity comparable to sound speeds, ~ 
a/ (o3/4p )[l + a,x/2az\ in the ferroelastic material at long 
times. An extremal textural profile £2(^7 ^) in Fourier space, 
if determined from an effective Lagrangian whose variation 
yields (3.36), will be shaped by this peaked structure. This 
suggests that textures will tilt diagonally and thereafter remain 
stationary and rigid; and that the most unstable or strongest- 
driven transient interfaces between phases will move/grow at 
a constant velocity close to the speed of sound, with bent 
domain-wall corners moving out, or by segment tips growing, 
along the diagonals. 

A fuller investigation of possible evolutions, with interme- 
diate states that can be sensitive to elastic and frictional pa- 
rameters, requires further work. Here we will only illustrate 
the rich variety of texturings, and later outline a tentative sce- 



FIG. 3: Square-rectangle (SR) case: grain-boundary motion under 
BG-type dynamics. The columns show (top to bottom) temporal 
sequences for time t of 10, 20, 42 and 70 'picoseconds (see Appendix 
B). The initial conditions are £2(r,t = Q),ez(r,t — 0) random 
around zero mean. Parameters (defined in Appendix B) are po = 
1, t = —0.25, and the material is 'hard', ai = 100, 03 = 210; 
with A' 2 = l,a' 3 = 0.1 = a'i . The time step is At = 10" 4 . 
Left column: The OP deviatoric strain £2 (f, t) under BG-equivalent 
dynamics (3.32), showing formation of twin-like regions separated 
by grain boundaries that are pushed out by domain-wall tip growth. 
Right column: non-OP shear strain ez(r,t), that is expelled from 
diagonal-domain regions, and concentrated in the pushed-out grain 
boundary regions (ei not shown). DiagnosticsEj at t — 70 were 
E = -1.08, (£2) = 0.0012, max-min e 2 = (1.33, -1.35); (e 3 ) = 
-7.9 x 10~ 6 , max-min e 3 = (0.062, -0.058). 

Figures 3 and 4 illustrate SR case simulations!! under the 
PG-equivalent dynamics of (3.32), with red/green/blue rep- 
resenting positive/zero/negative strains, with relative color in- 
tensities. Parameters are in the captions and quenches are into 




FIG. 4: Square-rectangle (SR) case: interface propagation under 
BG type dynamics. The columns show (top to bottom) temporal se- 
quences for time t in picoseconds of 40, 80, 160 and 1000. The initial 
conditions are £2 (r, t = 0) , e3 (r , t = 0) random around zero mean. 
Parameters (defined in Appendix B) are po = 1, r = —0.25, and the 
material is 'soft', a± = 10, 03 = 21; with A' 2 = 1 = a'3, a'i = 0. 
The time step is At = 0.002. Left column: The OP deviatoric 
strain £2(7% t) under BG-dynamics (3.32), showing domain walls 
propagating under the repulsive long-range compatibility potential, 
giving the impression of a 'zoom-lens' moving in. Right column: 
non-OP shear strain es(f, t), propagating outwards wjth interfaces, 
concentrated at corners, (ei not shown.) DiagnosticsEj at t = 1000 
were E = -1.4, (e 2 ) = 0.0044, max-min e 2 = (1.23,-1.23); 
(e 3 ) = -0.0168, max-min e 3 = (0.33, -0.49) . 



the martensitic regime r = —0.25. Figure 3 shows time se- 
quences of BG dynamics for the OP deviatoric strain and 
the non-OP shear strain. Note the grain-boundary-like re- 
gions that rapidly anneal out by tip growth, carrying the shear 
strain that is already expelled from the diagonal-domain re- 
gions. For these parameters, the propagation time t p < to, 
the unit-cell relaxation time defined later. Figure 4 shows 
the time sequences, for smaller elastic constants (and hence 
smaller shear mode velocities/larger propagation times), so 
t p > to- We now have domain walls propagating away from 
each other, giving the effect of a 'zoom-lens' moving in. Note 
the non-OP shear fronts moving with the OP walls (compare 
with Fig. 6 below). Such microstrupJure has been seen in 
FePd using phase contrast micppscopyo, and twinning waves 
have been found in ID modelsEl. 

For parameters as in Fig. 4 and for an intermediate temper- 
ature, t — 0.85, Figure 5 shows dynamical stress responses 
to applied static deviatoric stress P2(r) = Pq/[1 + (f/Vo) 2 ], 



FIG. 5: SR case: strain evolution under BG-type dynamics, 
with a fixed, time-independent, Lorentzian-profile local stress. 
The sequence (top to bottom), for time in picoseconds of t = 
40, 60, 76, 106 with the same parameters as Fig. 3, but now r = 
+0.85. Left column: Dynamic texturing of deviatoric strain £2(7', t). 
The system reduces the energy from the imposed single-sign strain 
by 'elastic photocopying', or adaptive screening of the long-range 
elastic interaction, generating higher multipoles that here propagate. 
Right column: The non-OP shear strain ez(r,t) follows the OP 
propagation. 



where rrj = 1 is the width and Po = 1 is a time-independent 
strength. The single-sign induced strain has a large energy, 
and the system elastically screens it by nucleating hierarchical 
opposite-sign elastic multipoles, with the propagation length 
setting a scale, and wave fronts moving out. A sinusoidally 
time-varying Po (t) can produce even more striking propagat- 
ing patterns. Similar 'elastic photocopying' was found previ- 
ously \md&c-TDGL dynamics, with the surface compatibility 
potentiaOa setting a domain-wall separation scalec2ci 

The last BG simulationS of Fig. 6 shows, for the TR case 
and (3.26), plots of the OP sheafcSap-with star-triangle pat- 
terns as found by other dynamicslijOLJ and seen experimen- 
tally in crystals of lead orthovanadate, Pfe^O^, that un- 
dergoes a trigonal to monoclinic transitioned. See also Fig. 8 
below. Once again, the e\ plot shows that equilibration in- 
volves expelling non-OP strains (except at OP corners, with 
ei global cancellation). 

Finally, we note that linearizing the BG dynamics (3.36) 
about equilibrium in the zero-damping limit yields the famil- 
iar wave equation. 'Textural phonon' spectra can emerge. 
For a finite Lq x Lq system we can also include surface- 
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FIG. 6: Triangular to centered rectangle (TR) case: nested strain tex- 
turing under BG dynamics of (3.26). The columns show temporal 
sequences, with time t in picoseconds of 10, 40, 250, 620. The initial 
conditions for the OP strains are £2(7% t = 0), £3(f*,t = 0) random 
around zero mean. Parameters (defined in text) are po = 1,t = 
-50, oi = 1000, a 3 = 2100; A' 2 = l,a' 3 = l.a'i = and the 
time step is At = 10 -3 . Left column: The OP shear strain £3(r, t) 
showing formation of nested star and triangle domains. Right col- 
umn: non-OP compression strain ei(r,t), concentrated near do- 
main corners, and expelled elsewhere. Diagnostics^ at t = 620 were 
E = -43.2, (e 2 > = -0.00316, (e 3 ) = -1.05 x 10~ 5 , max-min 
£2 = (4.45, -2.62), max-min e 3 = (3.63, -3.65). 



compatibility restoring forces, with a kerneH U surface ~ 
(a,3/aiLo)/\k\. The long wavelength OP strain oscillations 
then have velocities u>(k)/k = v(k) that are obtained by solv- 
ing 

v 2 ~ J- [< F" > +ai U c + a x u sur f ace ] , (3.42) 

where < F" > is a free energy curvature averaged with a 
probability distribution peaked at the equilibrium structure. 
For infinite systems the long wavelength spectrum is linear. 
For finite L , very long wavelengths probe the system sizeEl, 
and a; ~ (|fc|/Lo) 3 - This is the "dyadon spectrumEj of waves 
in twin-bands of martensite. Anomalies have indeed been 
observed^ in some ferroelastic phonon spectra. However, we 
do not pursue this conjectured explanation here. 



C. External stress and noise terms 

We now consider stress {p s (r)} and delta-correlated noise 
{ds(t)}> tna t modify the free energy as 

F -> F - ^2{pe£e + fjtst), 



/->■/- ^2(Piei + 9iCt), (3.43a) 

where noise correlations are 

< 9t(f> t') >= 2A' e fd er 5 Pr -,6(t - t'); 

< gi(r,t)gi'(f^,t') >= 2a! i T5 itV 8 f p5{t-t l ), (3.436) 

i.e. the bare noise is Markovian. Here T = ksT/Eo, and 
Eo is an OP elastic energy of Appendix B. For small stress 
and noise, a simple use of the substitution (1.4), justified by a 
detailed analysis, yields effective OP stresses and noises, 

Pi ot =Pt+ Pi c \ fff =9i+ 9t, (3.44a) 

where 

p/ e = & = J>S«, (3.44.6) 

i i 

with constants Su as in (3.24) and (3.34), respectively for the 
TR and SR cases. The total correlations are 

< 9t*(k,t)g*% t (k',i!) >= 2fA' v ,(k)6 Lg S(t - t'), 

(3.44c) 

where as in (3.28d) 

A' t ,e (k) = Stt'A'e + ^ a'^uS^ . (3.44d) 

i 

The elimination of non-OP strains thus induces cross- 
couplings, so non-OP stresses induce OP variations; noise 
correlations become spatially nonlocal; and different OP's 
acquire cross-correlated noises. The BG deterministic dy- 
namics then becomes the BG Langevin dynamics of (1.3). 
For the TR case the noise is delta-correlated in time. For the 
SR case with OP only, there is frequency-dependence, but 
this can be circumvented by considering {e2, e$} of (3.32) as 
our system. Thus in both cases we have two variables with 
Markovian noises. 



D. Fokker-Planck Description 

Langevin dynamics with delta-correlated Markovian noise 
can beCJ equivalently written in a Fokker-Planck (FP) de- 
scription. The set of 4L 2 o random variables labelled by 
a = {£, k} is taken to be 
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{x a } = {e 2 (k, t),s 3 (k, t)}; {v a } = {i 2 {k, t), i 3 (k, t)}. 

(3.45) 

for the TR case, with £3 — > e 3 for the SR case. The Langevin 
equations are 

x a (t) = v a (t), 

v a {t)=D^ a +f a {t), (3.46a) 
where the frictional force plus internal stress, or 'drift' term is 

1 d(F + F c 



DM, 



M(k) dx a * 



J2D {2) a ,a'V a ,, (3.466) 



where M(k) = 4po/k 2 is a strain mass. The Langevin noise 
correlation is 

< £ a (t)f a >(t') >= 2D^ a , a 'S(t - t'), (3.46c) 
with 'diffusion coefficient' 

D {2) a , a ' = 5 Z gTA' i4 ,{k)/M{k) 2 . (3.46d) 

The FP equation in Kramers form@ for the time- 
dependent probability P({x a }, {v a }, t) is 



dP/dt = LP= [L {1) + L<?>]P, 



(3.47a) 



where the Fokker-Planck operator for ferroelastics is a sum of 
drift and diffusion terms, respectively, given by 



-(2, 



E 



dx a *dx a i* 



D {2} , 

' y a. a' ■ 



(3.476) 



A formal solution for the probability in terms of the initial 
distribution isL3 



P{{x a } 7 {v a }, 1) = e tL P({x a }), {v a }, 0). (3.48) 

The FP operator carries the nonlinearity, symmetries, 
anisotropies, and long-range spatial correlations. Its eigen- 
values and eigenfunctions can be used to describe dynamic 
correlations. Since 'potential conditions' holco, the asymp- 
totic probability is a Boltzmann distribution, that (as can be 
checked by substitution) is P a = e -[( f + fC )+5L M ». 2 ]/ ? 
The free energy minima thus correspond to probability peaks 
in OP function space, highest for zero strain rate. For uniform 
OP this implies a triple-well free energy, but for nonuniform 
textures there will be a more complex free energy landscape. 
The multiple extrema are the TDGL asymptotic states. 

The FP formalism is convenient for discussing textural 
dynamics, metastability and glassy behavior; e.g. to deter- 
mine strain correlations thai—correspond to experimentally 
probed response functionsLfc3, or to calculate temperatum- 
dependent transformation rates through first-passage timesO. 
Both strains and strain-rates appear naturally,- as in phe- 
nomenological models of elasticity and plasticitycZl, that could 
thereby be given a microscopic basis. 



IV. BG DYNAMICS AS AN INHOMOGENEOUS ARRAY 
OF DAMPED AND COUPLED NONLINEAR OSCILLATORS 

Here, we consider a mechanical analog of nonidentical 
damped oscillators, that suggests a physical scenario for nu- 
cleation and growth after a temperature quench. We first re- 
view well-known damped oscillator results, to fix notation and 
terminology, and regimes of validity. 

A particle of mass M and friction parameter A, driven by a 
force —F' = —dF/dx = —Ax of spring constant A = F", 
obeys the underdamped equation 



Mx + A'x = -F', (4.1a) 
that can be written in an equivalent convenient form 



x = -(l/M)[dF/dx + A'x] 



(4.16) 



and the general solution is oscillations of exponentially de- 
creasing amplitude. With natural frequency ujq 2 = F" /M 
and relaxation rate r _1 = A' /M, the complex frequency 



^V^o-fo)- 8 "*^" 1 - (4-lc) 

Thus there is exponential decay (without even one complete 
oscillation) in the overdamped parameter limit of 



tqujq < 1, (4.2) 
when the inertial term is small at all times, and 



x(t) w -X[dF/dx{t)}, 



(4.3a) 



with A = 1/A', describes the damped behavior. Outside this 
regime it approximately describes the exponentially decaying 
envelope of the oscillations, at low frequencies/long times: 



w « A'/M; t » 2irM/A'. 



(4.36) 



Note that from (4. lc), and (4.3b), larger frequency/lower mass 
oscillations are damped out earlier. 

We now turn to BG-type evolution equations, that cannot 
be obtained by adding simple local inertial terms to the stan- 
dard overdamped dynamics like Model A or B dynamics for-a 
non-conserved or conserved order parameter, respectivelyCJ. 
However, defining an inverse mass 1/M(k) = 1/ p(k) = 
fc 2 /4po ~ k 2 the apparently unusual BG structure of (3.26) 
for the TR case, say, can be written as 

ii(k,t) = -(l/M{k))[d{F + F c )/det*{k,t) + A' e e e (k,t)], 

(4.4a) 

where for simplicity, a[ = 0. On comparing with (4.1b), the 
dynamics has an intuitively appealing interpretation. It is the 
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dynamics of a set of (nonlinear, coupled) oscillators, of (two- 
component) spring extension!^ ££(fc, t), labeled by k, with dif- 
ferent masses M(k) ~ 1/k 2 that are strongly dependent on 
the oscillator label k: heaviest near the origin, and lightest 
at the Brillouin-zone (BZ) edge. Similarly, the damping rate 
To(fc) -1 = A' /M(k) ~ k 2 is smaller for the larger masses. 
The spring coupling U c (k) acts equally over all the BZ oscil- 
lators, in a given direction k. Note that there is an intrinsically 
large range of damping times To ~ fc~ 2 , over orders of magni- 
tude. (This is reminicent of decay times at criticality at a sec- 
ond order critical point). From (4.4a) the k — ► infinite-mass 
oscillator is_a special case, with its initial velocity e(k — ► 0,0) 
unchanged!!]! 

Linearizing around equilibrium extensions, with effective 
curvature F" = F " + K k 2 = A(k), and with w 2 = 
A(k)/M(k), the complex frequency is 

J = V /c (fc) 2 -(r (fc))-2_ l I To (fc)- 1 . (4.46) 

The top panel of Fig. 7 shows that the angularly av- 
eraged strain structure factor S(k,t) =< |£T2(fc,i)| 2 >Sk 
or squared oscillator extension averaged over a shell 25k, 
is underdamped for small- fc oscillators, but is overdamped 
for larger-fc oscillators. The bottom panel shows asymptotic 
agreement between BG and TDGL dynamics. (See also 
Figs. 6, 8). The initial strains are nonzero only in a 5 x 5 
square region, and the oscillating strains show up as oscillat- 
ing colors (not shown). The bottom panel shows that for the 
SR case the moment < k 2 > takes on asymptotically the 
same value in the BG and TDGL dynamics. 

From (4.4b), the long wavelength modes are of necessity 
underdamped, with ujoTq ~ 1/k >> 1, a point made by Reid 
and GoodingQ. We will for simplicity consider the regime 
where none of the oscillators are in the overdamped regime 
(4.2), i.e. To(fc)wo(fc) >> 1 for all k, even those at the BZ 
corner k = y/2tv. The oscillating, damped texture has to set- 
tle down to some time-independent state, that from Sec. III.D 
will be a peak of the probability distribution, where the free 
energy derivative is zero. Then an 'envelope dynamics' as in 
(4.3a) for the OP strain is 



et{k,t) 



-A 



d{F + F c ) 
de e *{k,t) ' 



(4.5a) 



where A = l/A' e . However, the different-mass/friction oscil- 
lators can reach this late-time behavior only at different times 
tfc. Thus for times and lengths introduced through u ~ 2n/t, 
k ~ ir/L, there is a wavevector-dependent time scale t > % 
or a time-dependent length scale L < Lo(t) for textures to 
achieve the late-time regime, namely 

u « A'/M(k); t > t fc = 2n/[(A'/'ipo)k 2 }, 



L<L D {t) = (t/t D ) 1 / 2 ; tn^iApo/nA') 1 / 2 , (4.56) 
where t d is the time for relaxation across a unit cell. 
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FIG. 7: Top: Angularly averaged square oscillator extension or strain 
structure factor S(k, t) (defined in the text) versus time in picosec- 
onds; for S^-case nonlinear BG-dynamics oscillators, labeled by k. 
The oscillators are at high temperatures r = 4 (where there is only a 
single-well free energy), with parameters as in Fig. 4, and an initial 
condition as in text. The set of inhomogeneous strain oscillators of 
different mass ~ fc -2 and damping ~ fc 2 , are averaged over a shell 
of thickness 28k in the Brillouin zone. The heavier, less damped os- 
cillators with smaller fc = 0.525 (blue) oscillate, while the lighter, 
strongly damped oscillators with fc = 2.65 (red) are overdamped. 
The shells are 28 = 0.05 and 0.1, respectively. Bottom: Comparison 
of < fc 2 >= J2 k fc2 N 2 /E fc N 2 for BG and TDGL dynamics 
for TR oscillators. The BG line (blue) asymptotically merges with 
the TDGL line (red). Inset shows early time damped oscillations in 
the BG case and overdamped monotonic decay for the TDGL case. 



Taking the idea of a late-time equilibration length Ljj(t) 
as more generally applicable to the nonlinear case, a ten- 
tative picture emerges. While damped harmonic oscillators 
have trivial and identical final states, an inhomogeneous set 
of nonlinear oscillators with long-range coupling can have 
fc-dependent inhomogeneous final states {e(k, t — * oo)} = 
{e(k)}. The oscillations begin around the average value (say, 
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zero) of the initial states {e(k,t = 0)},for times t > tr> after a 
temperature quench. Because of the fc-dependent inhomoge- 
neous damping of (4.4a), the lightest masses near the BZ cor- 
ner, with labels tt > k > ir/Lrj(t) = fcu(i) will be the first to 
feel the final state attractor and begin oscillations around the 
final rather than the initial state. As time proceeds, the bound- 
ary shifts, and the circle of k < kp (t) 'initial-state' oscillators 
shrinks, while the number of \/2n > k > kjj(t) 'final-state' 
oscillators increases. Finally, all except the very smallest- k 
oscillators are on an equilibration path, and the larger-fc ones 
(corresponding to structures of time-dependent size Lo{t)) 
have already reached it. There is thus sequential-scale equili- 
bration, from the edge of the BZ inwards. For the SR caseEJ, 
the additional L p (t) propagation length enriches the scenario. 
The tentative scenario is consistent with what we have seen in 
our simulations for given parameters, although the nonlinear 
mode-coupling could induce more complex relaxation path- 
ways, in more general parameter regimes. 

Thus the unusual BG dynamics implies unusual elastic 
properties. Since p(k) ~ 1/fc 2 , long wavelength strains are 
kinematically blocked from decaying too early, and the scale- 
dependent damping and equilibration process starts at small- 
scale textures oriented by compatibility potentials, and then 
spreads to larger length scales, with associated non-OP strain 
expulsion. Consequently, materials classes governed by BG 
dynamics can have rich spatial patterning, metastability and 
glassiness, hierarchical multiscale microstructure, and a com- 
plex nonequilibrium elastic response. 



with the determinant (actually isotropic for A' 2 = A' 3 ) 



W 



^2^3 



^3 A 2 



(5.2c) 



Since diagonal elements, and the determinant, of the A ma- 
trix are nonzero, (5.2) implies that asymptotic textures ee(k) 
are determined by the extrema of the total free energy, that 
locate stable and metastable textural minima, as in the FP 
discussion of Sec. III.D: 



d{F + F c )/de\{k) = 0. 



In coordinate space, 



d(F + F c ) 



(5.3) 



(5.4) 



This is a time-dependent Ginzburg-Landau equation, with On- 
sager coefficients that are anisotropic and spatially nonlocal. 
For negligible non-OP friction parameter (a! x l A 2 — > 0), the 
Onsager coefficient matrix becomes both spatially indepen- 
dent and diagonal in OP labels, yielding an 'ordinary' TDGL 
equation with a constant, isotropic, and uniform friction: 



e e (f,t) = - 



1 d{F + F c ) 
~K t de e {r,t) ' 



(5.5) 



but still anisotropic and nonlocal in the elastic forces. 



V. LATE TIME/SMALL SCALE LIMITS OF DYNAMICS 

We now show (for all frictions a' it A' e nonzero) the 
underdamped TR and SR equations are approximated by 
TPGL-type equations. The validity o£ TDGL dynamics 
in any regime was recently questioned!!!] and an alternative 
overdamped dynamics, obtained by dropping displacement- 
acceleration terms in (2.8) was proposecOliZl. Appendix C 
shows they are equivalent. 

77? case dynamics: 

For the TR case, the late-time equations, from dropping 
OP strain-acceleration terms in (3.26), are with (3.28c), 



d(F + F c ) 
de* e (k,t) 



(5.1) 



This truncation is valid only for textures in a low fre- 
quency/large wave vector regime as in (4.5b). For a[ = 
(4.5a) follows, while for the general case, an inversion yields 



d{F + F c 



de*,(k,t) 

where the Onsager coefficient matrix is, with A' 2 = A' 3 



(5.2a) 



A's + a^Qh -aiQ 2 ,i<93,i 



/2,lQ3,l A' 2 



(5.26) 



SR case dynamics: 

For the SR case, from dropping OP inertial terms in (3.36) 
in fc, uj space, 

F 2 -iLu[A' 2 +r ] c (k 1 uj 2 , v 2 )]e 2 + U c (k, ui 2 , v 2 )e 2 = 0. (5.6a) 

This can be written as a generalized TDGL-type equation 
with retarded Onsager kernels A(r — r',t — t') in coordinate 
space 



dt'Y,^-^,t-t') Y^ F 2 - ^ 

de2\f ,v) 



Since the general TDGL structure is a 'current' proportional 
to a 'force', the Onsager coefficient X(k,ui) = l/[A' 2 + 
rf{k,u) 2 ,v 2 )] is like a dynamic 'conductivity'. We consider 
however, the low frequency ui — > 0, or asymptotic long-time 
limit, keeping the linear term in frequency, but neglecting 
the quadratic and higher order uj 2 , (uj/k) 2 frequency depen- 
dence in the kernels. The regime of validity is considered 
below. This yields an instantaneous nonlocal TDGL equa- 
tion as in the TR case, (5.4), with a single I = 2 order pa- 
rameter. In Fourier space, F c = \ u U c Q(k)\e 2 (k,t)\ 2 , 
where the retar ded compatibility potential reduces to the static 



expressionEcj of (3.41) and the Onsager coefficient \{k) is 
given by A(fc) = l/[A' 2 +77 c (fc, 0,0)]. 

For non-OP friction a[ 3 /a 13 — > 0, rf vanishes, and we 
again recover an ordinary (local, instantaneous) TDGL equa- 
tion as in (5.5) and I = 2 with compatibility forces remaining 
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nonlocal. Thus the model used in previous workup is a specific 
limit of the exact dynamics. 

In a regime similar to (4.5b) (with d' 2 = 
M ax[(a[ 03/0103), (a'j/ai) 2 ] ), the non-OP inertial de- 
lay (pq / a^){ui / k) 2 and frictional retardation ~ (cud 1 ) 2 
can be respectively neglected, yielding (5.5), for lengths 
L < Min[L D (t),L p (t)}, and times t > t f = {2-Kd') 1 / 2 . 
These simple heuristic estimates may not, of course, be 
strictly quantitative, but capture the diffusive aspect of the 
late-time relaxation. Figure 8, left column, shows that the 




the free energy density difference is then ~ 1/Lq. (A surface 
compatibility potential sets a domain-wall separation length 
scaleEl and would raise the degeneracy, favoring equal-width 
'true' twins.) The barrier between states with differing num- 
bers of walls is the cost ~ A\Lq of a fractional-length kink in 

the domain wall and the barrier crossing time ~ e t . 

Thus, in general the ALR potential can produce a multiple- 
minima free energy landscape, with nearly degenerate , 
differently-textured states separated by large barriers. Initial 
conditions or intermediate-state dynamical scales can lock the 
system into one of the metastable states. The possibility of 
multiple minima is consistent with recent analyses of models 
of competing short and long range interactionsEj. 

In general, the inclusion of noise in PG-Langevin simula- 
tions will enable the system to more easily find low-energy 
minima. We now turn to other symmetries, and give the dy- 
namics for the four crystal systems of 2D ferroelastic transi- 
tions. 



VI. OP STRAIN DYNAMICS FOR ALL 2D-SYMMETRY 
TRANSITIONS 



FIG. 8: Ordinary TDGL dynamics and OP textures: Left column: 
evolution for TR case with time t = 1,20, 150 and parameters 
ai = 1000, r = -50, A 2 = A' 3 = 1. Right column: SR case, 
ai = 100, 03 = 210, r = —0.25, A' 2 = 1 showing states evolved 
from different initial conditions, but the same parameters, illustrating 
nearly-degenerate multiple free energy minima discussed in the text. 

TDGL dynamicsEi for the TR case yields textures similar to 
the (longer-run) BG dynamics of Fig. 6. 

We note that the free energy can have several metastable 
minima with different microstructure, but closeby free en- 
ergy densities. Thus the nested stars of Fig. 6, under perturba- 
tion, yield rhombhohedral structures of slightly lower energy 
densityQ (that are also obtained directly with a different ini- 
tial random-number seed). The multiple-minima picture also 
emerges in TDGL simulations for the SR case of Fig. 8 where 
different random-number seeds produce three SR case quasi- 
twin textups, with diagonal-domains of different number and 
separatiorO. Such quasi-twins were previously obtained in a 
displacement picture by Monte Carlo simulatiorO The total 
free energies in each case are extensive ~ L$ 2 , while their en- 
ergy differences behave as the length ~ KqLq of the diagonal 
domain walls, (with compatibility cost from U c vanishing): 



Since the derivations above involved the generic strain mass 
tensor and the harmonic non-OP strains and did not involve 
the details of the anharmonic OP free energy, the BG dynam- 
ics structure will be the same for all types of 2D (and indeed, 
3D) transitions. However, the nature of the non-OP strains 
determines the symmetry of the compatibility kernels. 

In 2D there are N = 3 strains (compressional, deviatoric 
and shear) and one compatibility equation (7V C = 1). The 
distinct ferroelastic transformations (Fig. 1), order parameters 
and all possible symmetry-allowed combinations of strains in 
the Landau free energy, were found by Hatch et alH using 
the computer program ISOTROPY. These transformations fall 
into two classes determined by the nature and number of order 
parameters. (A) For either a two-component OP or two one- 
component OP's (N op = 2, n = 1), we have the TR case 
and the square to oblique (SO) case, respectively. (B) For a 
one component OP, (N op = 1, n = 2), we have the SR case 
with deviatoric OP; the square to centered rectangle (SC) 
case with shear OP; and the rectangle to oblique (RO) case 
with shear OP. In group theoretical symmetry notation, the 
transformations in Fig.l are: (a) P4mm to P2mm, (b) P2mm 
to P2, (c) P4mm to C2mm, (d) P6mm to C2mm, (e) P6mm to 
P2, and (f) P4mm to P2. The OP compatibility kernels are 
frequency independent (dependent) for case A (B). 



A. One non-OP strain, two OP strains (N op = 2, n = 1): 

• TR case, driven by combined deviatoric and shear 
strains e 2 , £3. This is studied in the text above. 

• SO case, driven by independent deviatoric e 2 and shear 
strains, £3. For the SO, there are two distinct OP's that 
drive the transition, and the OP free energy istj 
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F = ^Y / A 2 e 2 2 + B 2 e i 2 + A 3 el + B 3 ei (6.1) 

T 

where the constants are merely illustrative. However, the 
harmonic non-OP energy and thus the form of the dynamics 
and the kernels, are identical in both cases. 



B. Two non-OP strains, one OP strain (N op = 1, n = 2): 



Here, a\ = ai - a x Qi, 2 ; a 2 = a 2 - a x Q 2 ,i- The dy- 
namics is as in (6.3), but with the substitution ai,2 — > 
ai,2, and a[ 2 — > a' 12 , where a[ = a[ — a' x Q\,2\ a' 2 = 
a' 2 — a'xQ2,x- The extra term in (6.3b) is X = 
{ax £3 + a'x^d)Q^, i 2 Qi, 2- This again shows that verti- 
cal/horizontal domains are favored, but now the relative 
magnitudes of the original non-OP elastic constants a± 
act like symmetry-breaking fields, with x axis orien- 
tation preferred for ax > 0, and y axis preferred for 
ax < 0. The predicted behavior can be verified in sim- 
ulations. 



• SR driven by deviatoric strain, e 2 . This is considered 
in detail in the text above. 

• SC driven by shear strain, £3 The OP free energy now 
is as in (3.38a) but is Po({£ 3 })- The non-OP harmonic 
energy, f = h ^2 r a i e i + a 2 e 2 > an ^ tne dissipation is 
Rtot = [iEr a i e 'i 2 + a 2 e '2 2 +^ 3 e'3 2 ]- Thus the deriva- 
tion carries over, with the interchange 2^3, and, in 
particular, the compatibility factors and symmetry con- 
stants interchanged, O2 Q3, c 2 *-* C3 in the dynam- 
ical kernel of (3.40). The u> — > limit, analagous to 
(3.41), is the static result. 



ai U c (k) = (6.2) 

° V [l + (a 1 /a 2 )Q 2 , 1 2 ] 1 ' 

This is zero for k x = or for k y = 0, so domain walls 
are either vertical or horizontal, as can be confirmed in 
simulations. The dynamics for e 2 , £3 are: 

c 2 2 fc 2 

Poe'2 = — [ Q2,l<93,l(oi£3 + 

+ (aiQ 2 ,i 2 + a 2 )e 2 + {a\Q 2 ^ 2 + a' 2 )e 2 ], (6.3a) 



P0S3 = - ^ j [ Q2,iQ3.i(aie2 + aie 2 ) 

+ (P 3 + aiO^es) + K + aiOL)e 3 +X], (6.36) 

where the extra term X = in (6.3b). 

• RO driven by shear strain, £3. The non-OP free energy 
is harmonic in the sum and. difference of the compres- 
sion and deviatoric strainsO: 

/ = 5 E>+(ei + e 2 ) 2 + a_( ei - e 2 ) 2 } 

= \ Z)r{ a i e i 2 + a 2e2 2 + 2axeie 2 }, where a\ = 
a + + a_ = a 2 and ax = a + — a_ . Similiarly R tot = 
\ Sr-[^3^i+ a 'i e ' l2 + a 2 e '2 2 +2aj s :eie 2 ]. The derivation 
carries through. The static limit compatibility kernel is: 

iiU c (k) = Q 3 f [ ai _ + axQ.1,2) ■ 

\[1 + (ai/a 2 )Q 2 .i J / 

(6.4) 



VII. SUMMARY AND DISCUSSION 

There are two themes in this paper: firstly, a derivation of 
ferroelastic evolution equations for all 2D symmetries, and 
secondly, a demonstration of a strain-based (rather than a 
displacement-based) description of elastic phase transitions. 

We have, for the first time, derived the D > 1 underdamped 
dynamics for ferroelastics in terms of the order-parameter 
(OP) strains {e^} alone, showing that the evolution equations 
are of a generalized Bales-Gooding form. The strain-based 
derivation yields a wavevector dependent strain mass ~ 1/fc 2 , 
thus large-scale strains have greater inertia. The structure is 
that of an OP strain acceleration, e'i, proportional to a Lapla- 
cian acting on the sum of an OP-only stress and an OP- 
only frictional force. The stress and friction are strain and 
strain-velocity derivatives, respectively, of the effective free 
energy F({et}) + F c ({e{\) and effective Rayleigh dissipa- 
tion R({ee}) + R c ({ii}). These contain, in addition to di- 
rect local OP contributions (P and R), additional anisotropic 
and long-range contributions (P c and R c ) that emerge from 
eliminating the non-OP strains using St. Venant's compati- 
bility conditions. The kernels are explicitly evaluated for all 
2D symmetries. There are also compatibility-induced noise 
contributions, and this PG-Langevin dynamics of (1.3) or 
(3.46) is a central result. A Fokker-Planck equation (3.47) 
is obtained. The BG dynamics can be regarded as nonlin- 
ear, nonlocally coupled oscillators labeled by k, with unequal 
masses ~ 1/fc 2 , and dampings ~ fc 2 . The textures are the 
set of final rest positions {eg(k,t — * 00)}, with large-fc os- 
cillators (small-scale strain textures) equilibrating first. The 
late-time envelope dynamics that guides the damped oscilla- 
tions to this equilibration is of the TDGL form. This analog 
suggested an appealing picture of sequential-scale evolution 
for post-quench nucleation and hierarchical growth, account- 
ing for nonuniform textures. 

We adopt the strain picture in simulations, with strains as 
the basic variables on sites of a reference lattice, driven by 
symmetry-allowed terms of the Ginzburg-Landau free energy, 
and by anisotropic, symmetry-specific, long-range compati- 
bility forces. The free energies are in a standardized form, 
with dimensionless parameters related to experiment. Simu- 
lations show that the BG dynamics has rich texturing prop- 
erties, including repulsive velocity-resonant compatibility po- 
tentials that can drive interfaces at (nearly) sound speeds. 
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We now place our results in perspective with some of the 
other models thai,have earlier provided valuable insights. 
Baus and-. Love ttt3 invoked the 19th century work of St. 



and. 

VenantEE3 on the compatibility condition for the strain tensor, 
in the context of surface tension in liquids. They considered 
strain as the basic variable in the argument, and noted it might 
be useful in elastic solids. We similiarly work | '" thp strain pic- 
ture, differing in this respect from previousLxja simulations 
that work with u gradients, i.e. in the displacement picture. 

In an interesting and important paper, Kartha et alo per- 
formed Monte Carlo simulations to find static textures. They 
used the SR free energy that is sixth power in the deviatoric 
strain order parameter and harmonic in the compressional and 
shear non-OP strains, and the simulations were in terms of 
the displacement, so effectively V — V({u}). Since strik- 
ing textures, like unequal-separation diagonal domain walls 
(and tweed) were obtained, they attempted to understand these 
u simulation textures by using compatibilityLxj to eliminate 
non-OP strains in the free energy, plotting (static) compati- 
bility potentials. However, these effective strain-strain corre- 
lations were not directly used in the simulations. Such ex- 
plicit implementation of the compatibility forces was done 
in a TDGL strain dynamics, where quasi-twins, the elastic 
Meissner effect (expulsion of non-OP strains) and tweed were 
investigated!!!! Other JXDGL displacement simulations-ii 
vestigated tweed alonetl (The tweed terms consideredLiHl 
were all different.) 

Our work has the same Lagrange-Rayleigh starting pointBEa 
as Refs. [l6|JT^ , that focus early on in their argument on an 
overdamped-displacement dynamics. We follow a different 
path and derive an OP-only underdamped strain dynamics, 
finding it is of a generalized BG form; and that the late-time 
limit is a TDGL-type equation. 

Our approach differs from an underdamped dynamic so 
for the SR case, that did not explicitly consider non-OP 
strains; and phenomenologically added a static anisotropic 
long-range potential between squares of strains to explain 
acoustic signalsE! We derive the BG structure from a La- 
grangian with both non-OP strains and compatibility con- 
straints, and find a retarded anisotropic long-range force in 
terms of the OP strains themselves (and not their squares). 
Our dynamics in the TDGL limit also differs from a TDGL 
dynamicsffil, where strains have been eliminated in favor of 
morphological profile variables, r] a (r, t), with a labeling the 
structural variants. The static potentials between squares of 
the morphologicaLvariables were obtained from elastic fields 
due to inclusionsEZl 

Our approachis in the spirit of the Landau description of 
phase transitionsEj: working with the order parameters as the 
basic and physically relevant variables, and focusing on the 
order parameter symmetries (as encoded in the compatibility 
factors), as the source of ferroelastic static and dynamic tex- 
turing. 

Further work could include a detailed understanding of 
2D nucleation, growth and interfacial profiles; extensions 
to 3D symmetries such as cubic to tetragonaltil (N op = 2, 
n = 4, N c = 6); generalizations to include defects, in a 
broader 'strain elastodynamics' framework; making contact 



with phenomenologies of plasticity; simulations and calcu- 
lations of experimentally measurable strain correlations and 
nonlinear susceptibilities: and exploring a hierarchical sce- 
nario for shape memoryEl 

The symmetry-specific, compatibility-focused under- 
damped ferroelastic dynamics for the strain order parameter 
encode, in their very structure, the possibility of an evolution- 
ary textural hierarchy in both space and time, and a tendency 
for interfaces to be driven at sound speeds, explaining some of 
the fascinating but puzzling features of martensitic dynamics. 
The dynamical equations could be applied to a wide variety 
of textural evolutions that include improper ferroelastics, 
leading to a deeper understanding of many materials of 
technological interest such as ferroelectrics, magnetoelastics, 
colossal magneto-resistance manganites, superconducting 
cuprates, and shape memory materials. 
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APPENDIX A: SQUARE TO RECTANGLE TRANSITION 
DYNAMICS 



We need to demonstrate explicitly that the same dynamics 
results whether displacements or strains are treated as the ba- 
sic independent variables. We derive, in two self-contained 
subsections the (same) underdamped dynamics for the SR 
case, (a) by varying the displacement, and (b) by varying the 
strains subject to compatibility. 



1. Variation of displacements 

The (n = 2) non-OP strains are here the compressional 
(ei) and shear strains (e^), and the N op = 1 deviatoric strain 
(e 2 ) is the OP: 

— = \ {A X U X + AyUy) , — = l (A X U X - AyUy) , 



e 3 



1 



(A X Uy + AyU X ) , 



c 3 2 

with the free energy V = f + F, where 



(Al) 



(A2a) 
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The Rayleigh dissipation function is 

R tot = 1^21^1 + a '^+ A '^- (A2b) 

T 

From (A2) and varying with respect to u as in (2.1) and (2.2) , 

Po il x = ^{A x [fi + Ri] + \[fa + Rs] + A X [G 2 + R 2 ]}, 

(A3a) 

Po u y = \{A v [h + + A x [f 3 + R 3 ] - A y [G 2 + R 2 }}. 

_ (ASb) 
This is the result of Ref. [n], where 



, _ df 8F 

J 1,3 = Cl,3 ~ 5 Lx2 = C 2 — — ; 

9ei 3 oe 2 



Defining = [a 3ul - p (^-) 2 ]/a lLJ and B u = 1 + 
buiQi/Qz) 2 , we obtain 



ei(fc,w) = (Qi/Q 3 )b w e 3 (A;,a;), 

e 3 (fc,w) = -{Q 2 /Q 3 B LO )e 2 (k,u). 
Using (A6), and (A7) 



(A7) 



P0 Lo 2 e 2 = ^-k 2 [F 2 - iujA' 2 e 2 + a lw (Q 2 /Q 3 ) 2 (^/ 'B u )e 2 ]. 

(A8) 

The dynamics is written out in BG form and discussed at the 
end of Sec. III. 



2. Variation in Strain 



dR tot dR tot 
Ri,3 = ci,3— — ; R 2 = c 2 - 



dei 



di 2 



The underdamped strain equations with A 2 = 

A 2 ,D 2 ee A 2 - A 2 , are: 



(A3c) 
A 2 + 



Po ei = j[A 2 (/i + Ri) + 2A x A y {h + R 3 ) 



By varying (3.1) with respect to ei j3 (f, t), e 2 (r, t), we ob- 
tain N + N c = 3+1 = 4 equations: 



P22^2 + P2iei 



5F 
"SeX 



Q 2 A-A' 2 e 2 , (A9a) 



= -a 3 e 3 - a 3 e 3 - Q 3 A, 



(A9b) 



+D 2 (G 2 + R 2 )}, 



(AAa) 



{puei + /?i2£2} = -aiei - a[ei - QiA - a' 3 e 3 , (A9c) 



Po£2 = j[+D 2 (.h + i?0 + A 2 (G 2 + R 2 )l (A4b) 



Po e 3 = ^pA.A.jC/i + i?i) + A 2 (/ 3 + i? 3 )]. (AAc) 
We also have the compatibility relation of (3.6) 

Qiei + Q 2 e 2 + Q 3 e 3 = 0, (A5) 

where Qi = A 2 /a,Q 2 = -D 2 /c 2 ,Q 3 = -2A x A y /c 3 . By 
taking derivatives of (A4) we see that (A5) is identically sat- 
isfied. Thus we can take (A4a),(A4c) and (A5) as the three 
equations to determine the three strains, since then two of the 
equations are linear. Fourier transforming these three equa- 
tions, we obtain for ei. 3 (fc, u>), e 2 (k, u>) 

Pqu?e 2 = i[cic 2 (fc 2 - ky)a lu] ex + c 2 2 k 2 (F 2 - iujA' 2 e 2 )}, 



(A6a) 



Pquj e 3 = -[ci^kxkyai^ei + c 3 k a 3uJ e 3 ], (A6b) 



Qid + Q 3 e 3 = -Q 2 e 2 . 



(A9d) 



In (k, us) space for ei i3 (fc, ui), e 2 (k, ui) and with ai„ ee a\ — 
iojai, a 3u ee a 3 — ioja 3 , there is one nonlinear equation, 

SF 

u 2 p 22 e 2 + to 2 p 2 iei = +— — - + Q 2 A - iujA 2 e 2 (A10) 
6e 2 

and n + N c = 3 linear equations 

= -a 3lJ e 3 - Q 3 A, (Alia) 



( f\ ) e 3 = ~ a l" e l + Qi A = °> ( Allb ) 



}i&i + Q 3 e 3 = -Q 2 e 2 , 



(Allc) 



where we have used a relation as in (3.21). From (Alia) we 
have the Lagrange multiplier 



A(fc,w) 



a 3 ^e 3 



(A12) 



Qiei = -Q 3 e 3 - Q 2 e 2 , 



(A6c) 



where F 2 = dF/de 2 (k,ui); a luJ ee ax — jwa^; a 3w ee a 3 — 
iwa 3 andQi = k 2 /c 2 ,Q 2 = (A: 2 - k 2 )/c 2 ,Q 3 = 2k x k y /c 3 . 



and hence the non-OP strains in terms of each other, ei 
+&w(§ 1 )e 3 . Finally in terms of the OP, 



e 3 = -(Q2/Q 3 )e 2 /5 a 



(A13a) 
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where 



bu = 
Then 

where 



[«3^ ~ pg{2uj/kf 



i+bUQi/Q 3 ) 2 - {Am) 



4k 2 



e 2 = + [F 2 - iu)A' 2 e 2 + We 2 ] 



(AUa) 



W 



2f,2 



('2 



''[P22 - 



>1W2 



Q 2 2 



P. 



(Q 2 /Q 3 ) 2 - 



P. 



(AUb) 



We use the definitions (3.12) 

Pll = (^/^^^(dQx/csQa) 2 ; 
(ip /k 2 c 3 2 )(Q 1 Q 2 /Q 3 2 )pi2 and (Qi 2 - 
(c 3 /c 2 ) 2 = 1/2 to obtain: 



of 



P22 



P21 

Q2 2 )/Q 3 2 



C 2 

P0 Lo 2 e 2 = -^-k 2 [F 2 ~ iuA' 2 e 2 + a lw (Q 2 /Q 3 ) 2 (^/P w )£ 2 ]. 

(Alb) 

Comparing with (A8), we see that the same dynamics results, 
whether displacement or strain is regarded as the independent 
variable, with compatibility enforced in the latter case. The 
BG evolution equation is discussed at the end of Sec. III. 

A useful intermediate equation, obtained by substituting e\ 
in terms of e 3 into (A10) and using (A12) for the Lagrange 
multiplier, expresses the equations as a coupled set for the OP 
and non-OP shear. This dynamics is given in (3.32), and is 
entirely equivalent to the OP-only retarded equation of (A15). 



APPENDIX B: SCALING OF FREE ENERGY 

Here we.-eeneralize the free energy scaling of Barsch and 
KrumhanslLI to scale the Lagrangian such that the Landau 
free energy is in a standard polynomial form; all the parame- 
ters are dimensionless; strains are of order unity and times are 
scaled in a characteristic time unit. 

The symmetric strain tensor is related to displacement 
derivatives through 



1 Mfl + dUJfU + y dUp dUp 

2 dx v dx^ g Xfi g Xv 



(PI) 



and we retain the 'geometric nonlinearity', showing when it is 
negligible later. The symmetry-adapted compressional (</>i), 
deviatoric (</> 2 ) and shear (0 3 ) strains are 



h/c 2 = —\ 



fo/ c 3 = ^(^xy + <t>yx), 



(P2) 



where c±, c 2 , c 3 are symmetry-specific constants. The free en- 
ergy depends on the order-parameter through a Landau term, 
and has non-OP and gradient contributions: 



F Bj^ anc [ au -\- F non -\- Bq tcl( i- 
The order parameter Landau energy is 



(P3) 



F Landau = f d\^^B {2 \T)4>}+F poly \, (54a) 

where B^ (T) = Bq (T— T c ) vanishes at a characteristic tem- 
perature T c , and the F po i y is a temperature independent OP 
polynomial with powers higher than quadratic, that is different 
for the SR and TR cases. The non-OP contribution is 

F non = \Jd 3 rY / B 
while the gradient term is 



(B4b)) 



F G rad = I d\Y,\ K ^M 2 - (B4c) 

We assume uniformity in the z-direction, of thickness h, 
and pass over to a 2D reference lattice of lattice constant a Q 
(square or triangular for SR and TR cases), thus unit-cell 
strains are tensorial variables on the dual lattice. Derivatives 
are converted to discrete lattice differences, and displacements 
are scaled in a Q . Thus 



_d_ 

dx,, 



1 



A M ; / d 3 r^ha 2 J2 



(P5) 



We further scale strains in a typical value A chosen for con- 
venience later, so a scaled strain tensor is then defined 
through 

M „ = AP M „ = ^ [{A M u y + A„u M } + A ^ A^Up A„u p ] , 

(P6) 

so V ' fJJ v — ► AApMy where u is a dimensionless scaled dis- 
placement. The scaled OP and non-OP strains are 



<f>t = \ee, (j>% = Ae 4 . 



(B7) 



The free energy (B3) can be written in terms of these scaled 
strains as 



F = F a ({s e }, t) + /({ej) + F grad ({Ae e }), (B8) 
where 

Fnon/E = f = - ^2 a i e h 



Farad/ Eq = Fgrad — ^ ~^{^ £ lf 



(B9a) 



(B9b) 
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Here we multiply and divide by an energy density Dq chosen 
later, to get the overall energy scale E Q = ha 2 D Q (this drops 
out in the end). The scaled parameters are 



5,A 2 KX 2 
a t = ho = 



Da 



The scaled OP free energy is 



KD ) 



(BIO) 



F Landau /E = F = ^(t - l)e 2 + F 00 , (511a) 



where Fqo = F po i y /En + J2^i£i 2 - The scaled temperature r 
contains a physical temperature Tq that fixes where r = 1: 



_ (t-t c ) ,5 a 2 i 

T= — — , ) = 



(T -T c ) ! v D ' (T -T c ) 



(5116) 



We now consider the Si? and Ti? symmetries separately. 
SR case scaling: 

For the SR case, iV p = 1, there is only the deviatoric 
strain, as <j>2 as the order parameter and the polynomial is 



F poly = I d 3 r{-B^4 + B^cj> 2 e: 
In scaled form, the term in (Bl la) is 



(£12) 



units of ergs /cm 3 , B\ 
Co = 1.7 x 10 12 , D' 



2 

ergs 



-- 1.4 x 10 12 , 5 2 = 2.8 x 10 12 , 
3 x 10 17 , 4- = 2.5 x 10 11 , and 



B = 2.4 x lQ J c J'l b gK - This gives, from (B6), the typical 
strain value A = 0.02; elastic constants a\ = 155 f» 03/2; 
OP variation scale y/Ko w 5; an elastic energy density 
= D = 3.8 x 10 6 |^|, and a temperature separation 
To — T c = 7 Kelvin. (The magnitude of Do corresponds to 
a magnetic energy density ^ for fields ~ 1 Tesla.) Note 
that A << 1, thus in (B6) we may drop the 'geometric 
nonlinearity' and work with the linear Cauchy strain tensor, 
as in the text, that satisfies the simple St. Venant compatibility 
condition. External stresses (eg. compressional) enter as 
F — > F + J2 r Pi e i w i tn scaled pressure pi = 1 correspond- 
ing to £a = 0.02G5a. 

Inertial and Damping terms: 

Using the same kind of transformations, the dimension- 
less inertial (T) and damping (R tot ) contributions to the La- 
grangian are: 



p u 



(518a) 



and 



5 00 = 5Z[£2 2 - C £ 2 4 + e 2 6 ] , (513) 

r 

where we have factored out an energy density D = B^X e 
so the coefficient of the sixth order term is unity, and Co = 
B ( 4 ) \ 4 /D . Now we choose A to fix Co so that for three (a = 
+, — , 0) roots e 2 a at r = 1, the conditions of degeneracy 

F 00 ({4 a) m = F ({4 a) (l)},r= 1) = 0, (BU) 
and normalization 

5>i ±} (i)] 2 = i (5i5) 

I 

are satisfied, and hence determine the typical strain A and all 
scaled parameters. For the SR case, 

C = 2; A = (B< 4 V2B (6) )2, D = (flW /2) 3 /(B (6) ) 2 , 

(516) 

with the Co 'completing the square' in F o and thus the OP 
free energy is 

F = J2(T-l)4 + el(e 2 2 -l) 2 . (517) 

r 

The t = 1 roots, = ±1, then manifestly satisfy the con- 
ditions (B14)and (B15). 

To get an idea of parameters, we use FePd shape memory 
alloy valuesEj, for the SR case with energy densities in 








= 2X> / i«i + °3*3 + ^' (5186) 

r 

where the dots are dimensionless time derivatives, and we in- 
troduce a characteristic time unit t a , with dimensionless den- 
sity po and friction coefficients a[ 3 and A' 2 defined as 

Pm X 2 al < 3 B' h3 B' 2 /B Q tl 
*o5 ai l3 5i j3 ^ (T Q - T c ) 

Here we have used (BIO), (Bllb) for a h3 and (T a - T c ). 
Wave propagation crossing a nanometer in a picosecond 
corresponds to a sound speed of 1000 m/sec. We take t a to 
be of the order of inverse phonon frequencies and the scaled 
friction coefficients a[ 3 , A' 2 to then be less than unity. With 
a mass density of p m = 10 gm/cm 3 , lattice constant a a ~ 3 
A and t ~ 10~ 12 sec the dimensionless density, or ratio of 
kinetic and elastic energy densities is po ~ 1. We will work 
with parameters po = 1, Kq = 1, to = 1 picosecond, a fixed 
ratio 03/01 = 2.1, and a\ = 100 or 10, with a[,A' 2 ,a' 3 as 
unity or less. 

TR case scaling: j-j— , 

For the TR case we start fromtll3 a free energy density in 

F po iy ~ -B^^ 3 - 30201) + 5< 4 ) (0 2 + 0l) 4 . We follow 
the same procedure as above : (i) scaling strains as in (B7); 
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(ii) pulling out a common factor Dq = B^X 4 ; (iii) choosing 
A in C = B^X 3 /D to satisfy (B15). This gives 

Co = 2, A = B^/2B^\ D = (B (3) /2) 4 /(B^) 3 . 

(B20) 

The scaled free energy is then 



Foo = ^(£2 2 +£3 2 )-2(^-3£2£i) + (£2 2 +£3 2 ) 2 (521a) 



(a) TR case displacement truncations: Dropping and 
keeping A m m„ in (2.8) yieldal3 equations that 
in Fourier space are 

a[k x ei + A 2 k x i 2 + A' 3 k y i 3 = -{aik x ei + k x F 2 + k y F 3 }, 

(C2a) 

a^kyei - A 2 k y i 2 + A 3 k x e 3 = -{a x k y ex - k y F 2 + k x F 3 }. 

(C26) 

Using compatibility (3.18d) we can eliminate 



= £ 3(1 + 2£ 2 )[£ 3 2 - ^(1 - £2 2 ) 2 ] + [£2 2 + £ 3 2 - I] 2 - 

(5216) 

Here the second form explicitly displays the (B14), (B15) con- 
ditions, for r = 1 , when the roots are (£2 i £3^) = (1, 0); 
(£2 < - ± \ £3^) = (— |, ±^), and lie on a unit circle. 

Connection to other scalings: 

The SR free energy of Ref. [l7| can be written as our 
(B9), (B17), (B18b) by a scaling of strains in a = 10 6 , 
giving an overall factor \o? absorbed in the TDGL time. 
(The elastic and frictional constants are then half our Oj, aj.) 
A similiar TR scaling in a = 10 3 yields our (B21), with 
a common factor of i relating elastic/friction constants^. 
Similar scaling can be performed for other symmetries. 



ei = -[Q 2 (fc)£2 + Q3(fc)£ 3 ]/Ci(fc), 



(C3) 



h 2 —1c 

can be written in matrix form as 



where Q2/Q1 = (^p^);Q 3 /Ci = nF*-. Further, (C2) 



M 



£3 J \k y -k x J \F 3 + F§ J ' 



F 2 + F< 



(C4) 



where Mu' is defined as M 22 = k x [A' 2 — a[Q 2 /Qi\; 
M 33 = -k x [A' 3 - (fcv/ferJaiCQa/Qi)]; M 23 = k y [A' 3 - 
oKfe/^JCQs/Oi)]; M32 = k y [A' 2 + ai(Q 2 /Qi)]. Here 
F 2 , F 3 are chosen to match the RHS terms of (C2), so 
that F$ = oiJ(Q|/Qi 2 )e 2 + (C&Qa/Qi 2 )^], and F 3 C = 
ai[(Q2C3/Qi 2 )£2 + ai(C3/Ci) 2 £3] and can be written as 
derivatives of a compatibility potential F c : 



F 2 ,3 + ^2,3 — 



d(F + F c ) 

d£2,3*(k) ' 



(C5) 



APPENDIX C: TRUNCATED DISPLACEMENT DYNAMICS 
AND TDGL EQUATIONS 

We have shown that the BG equations, dropping strain- 
accelerations in an uj — k regime, yield TDGL equations. 
On the other hand, dropping displacement-acazlex&tions in 



(2.8) or A3 yield equations stated to be differentOlJ from the 
TDGL form. In this Appendix, we demonstrate their equiva- 
lence to TDGL. 

More generally, the truncation is like dropping all the iner- 
tial terms in the Lagrange-Rayleigh equations (3.1), (3.2): 



(CI) 



where F c = ^a 1 J2kU c u'£e(k,t)e*i'(k,t) with 
U c u > = QiQe>/Q 2 i- Inverting the matrix M of (C4) 
yields ,(5.2), and then the TDGL link is as in the text. For a 
choicdl3 a[ = 0, one gets the ordinary (local, instantaneous) 
TDGL equation, namely (5.5). 

SR case displacement truncations: 

Similiarly for the SR case, dropping displacement acceler- 
ations compared to gradients of displacement velocitiestZl in 
(A3) yields the Fourier space equation, with a iuJ = a\ — iuoa 1 ^ 



( — )oi w ei + ( — )a3o;(^-)e3 
c 2 c 2 k x 



-iu!A' 2 + F 2 ], (C7a) 



;eE 



given as a balanceM between derivatives of the stress ten- 
sor o^y = SL/SEnj, and the damping force tensor a' = 
SR^SE^. ClearlyO, one cannot proceed in 2D by simply 
dropping the d v derivatives in (CI). Such a procedure would 
give a^y = —cr', which is not quite correct. For the SR 
transformation, for example, this is e\ t3 = —{a-y^/a'y 3)61,3; 

and £9 = — tttst-, where F ~ £9 is just the local triple-well 

A A' 2 de 2 ' 

free energy. This is the kind of overdamped dynamic sQ, with- 
out compatibility contributions, that would emerge if all the 
strains were independent. 

We now show that the truncated displacement dynamics 
(CI) is in fact a TDGL dynamics. 



(-)o lw e x + (-)a 3 ^)e 3 = [-iuA 2 + F 2 \. (C76) 
c 2 c 2 k y 



Then compatibility ei + (gf-)£2 + (^) e 3 ~ gives 



Q 3 



W3 0-lu> 



e-3 



'! + (^)(§) 2 



]£2- 



Hence 



(£78) 
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FIG. 9: Evolution from single-site initial condition for SR case un- 
der ordinary TDGL dynamics. Times are t = 0.13, 0.16, 0.18, 0.24 
picosecs, with At — 10~ 4 . The initial condition is £2(r,t — 
0) = 0.0001 at a single site and zero elsewhere. Non-OP fric- 
tion constants are a[ — = a' z . Left column: 'Soft' case, 
ai = 2, as — 2,r = — 50, A 2 = 2. Right column: 'Hard' case, 
ai = 2000, a 3 = 2000, r = -50, A' 2 = 2. (Since colors are 
relative, the background changes, with changes in evolving average 
intensity.) 



-iwA 2 e 2 = -F 2 + 



£2- 



(C9) 



Separating real and imaginary parts of the square brackets, 
and using the notation of (3.35), 



-A(fc,w)[F 2 + ai(7 c (fc,u; 2 ,0)£ 2 ] 



(C10)) 



where now X(k,u>) = l/[A 2 + a[r] c (k,uj , 0)]. This is man- 
ifestly a generalized TDGL as in (5.6a), but with a partial 
truncation of the kernel, that drops the (resonant) inertial delay 
terms (po/ '03) (a; /k) 2 and keeps only frictional retardation. It 
is thus valid for a narrowly restricted, intermediate length and 
time regime roughly estimated as Lr)(t) > L p (t) > L and 
t < tf (see Sec.V). Taking a more well defined limit of late 
times (dropping all oj 2 ) the static kernel U c (k, 0, 0) becomes 
a good approximation, and the nonlocal TDGL (5.4) with. 
1,1' = 2 holds. For a[ 3/01,3 << A' 2 as for 'hard' systemstll 



this collapses to a local TJJ.GL, as in (5.8). Thus overdamped 
displacement equations 



are TDGL equations in disguise 
and not a new dynamics. As expecteiL-simulation of TDGL 
equations produce the same texture soL3 from similar initial 
conditions. 

Figure ^Lshows ordinary TDGL simulations with an initial 
conditionO of strain nonzero at a single point and zero else- 
where. This is like a stress applied at a single point and then 
removed. Flower-like or diagonal-cross textures similiar to 
Ref.[l7| are obtained, for both 'soft' and 'hard' materials. 

Reference |l7| could not reproduce a TDGL structure given 
in Fig. 4b of Ref. [To], part of a multi-panel figure that dis- 
plays stress effects for soft materials. In the intermediate- 
temperature (1 > t > 0) phase, Ref. [To] considered the effect 
of two Lorentzian deviatoric stresses, fixed and continuously 
maintained. The motivation was to see if a stress seed analo- 
gous to a seed crystal in a supercooled melt, could yield stress- 
induced martensitic twins even for positive r = 0.3. The Ref. 
17 simulations did not have the applied constant stress of Ref. 
1C, and moreover, were at very low temperatures t = —50. 
We conclude that the difference in results is due to a differ- 
ence in states investigated, and not a difference in dynamics, 
which is TDGL-like in both cases. 
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